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1.  Introduction 

The  possibility  of  the  expansion  of  clusters  of  galaxies  has  been  considered 
from  two  points  of  view.  First,  if  one  assumes  that  the  redshift  of  galaxies  reflects 
velocity  of  recession,  the  question  arises:  do  the  members  of  clusters  of  galaxies 
participate  in  the  general  expansion  indicated  by  the  law  of  redshifts,  with  the 
clusters  themselves  expanding?  In  a  study  of  this  problem,  Neyman  and  Scott  [1] 
developed  a  statistical  kinematic  test  of  the  stability  of  systems  of  galaxies.  In 
1953  the  test  was  applied  to  the  twelve  galaxies  in  the  Coma  Cluster  for  which 
the  necessary  data  (radial  velocities  and  magnitudes)  were  available.  There  was 
no  indication  of  expansion. 

More  recently,  the  possible  expansion  of  systems  of  galaxies  was  considered 
from  a  different  point  of  view.  The  suggestion,  emphasized  particularly  by 
Ambartsumian  [2],  is  that  some  systems  of  galaxies  are  disintegrating  rapidly. 
This  suggestion  is  motivated  partly  by  the  observation  that  the  average  mass 
of  galaxies  derived  through  the  application  of  the  virial  theorem,  valid  for  stable 
clusters,  sometimes  appears  to  be  considerably  larger  than  the  estimates  obtained 
by  other  methods,  for  example,  by  the  study  of  double  galaxies  by  Holmberg  [3] 
and  by  Page  [4].  The  discrepancy  has  often  been  interpreted  as  indicating  the 
presence  of  large  masses  of  intergalactic  material.  However,  some  investigators 
consider  that,  at  least  in  some  clusters,  the  discrepancy  is  too  large  to  be  ex¬ 
plained  entirely  by  dark  matter,  although  it  may  indeed  be  present;  they  con¬ 
clude  that  these  clusters  are  unstable  and  the  virial  theorem  is  not  applicable. 

The  basic  assumption  underlying  this  argument  is  that  the  various  galaxies 
for  which  divergent  estimates  have  been  obtained  by  different  methods  have 
average  mass  of  the  same  order  of  magnitude.  This  may  well  be  true.  However, 
a  priori  it  is  not  evident  that,  for  example,  the  members  of  small  systems  of 
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galaxies  studied  by  Page  must  be  of  the  same  mass  as  the  galaxies  of  the  giant 
clusters.  In  the  studies  concerned  with  the  estimation  of  masses  of  galaxies  and, 
in  particular,  with  the  inconsistency  of  estimates  obtained  through  the  applica¬ 
tion  of  the  virial  theorem  [5]— [11],  a  considerable  number  of  other  assumptions 
are  used.  These  assumptions  vary  from  one  study  to  another  and  occasionally 
lead  to  diametrically  opposite  conclusions. 

In  these  circumstances,  the  possibility  of  approaching  the  problem  of  stability 
of  systems  of  galaxies  from  an  entirely  different  point  of  view,  with  the  test 
criterion  published  in  1954,  totally  independent  of  any  assumption  regarding 
masses  of  galaxies  and,  indeed,  independent  of  dynamical  considerations  of  any 
kind,  presents  some  interest.  This  interest  is  enhanced  by  the  fact  that,  in  con¬ 
nection  with  a  contemplated  joint  study  of  the  Coma  Cluster,  Mayall  [12] 
undertook  an  extensive  observational  program  yielding  the  radial  velocities  and 
the  apparent  magnitudes  of  a  number  of  presumed  members  of  the  Coma  Cluster. 
The  number  of  objects  for  which  both  these  data  are  available  has  increased 
from  12  in  1953  to  50  galaxies.  Also,  combining  the  observations  of  various 
investigators,  currently  we  have  usable  data  for  a  number  of  other  clusters  and 
groups.  It  is  true  that,  with  the  exception  of  the  Virgo  Cluster,  in  each  case  the 
number  of  available  galaxies  is  much  smaller  than  for  the  Coma  Cluster.  How¬ 
ever,  there  is  the  possibility  of  combining  the  data  to  perform  a  summary  test. 

Actually,  two  summary  tests  are  performed.  Both  test  the  hypothesis  that  all 
the  clusters  and  groups  considered  are  stable.  However,  the  first  test  criterion 
was  deduced  to  be  particularly  powerful  against  the  alternative  hypothesis  that 
all  the  systems  of  galaxies  studied  are  expanding  at  the  same  rate.  The  second 
criterion  tests  the  hypothesis  of  stability  against  the  alternative  that  at  least 
some  of  the  systems  studied  are  unstable:  some  may  be  expanding  and  some 
contracting  not  necessarily  at  the  same  rate.  The  mathematical  assumptions 
concerning  the  properties  of  systems  of  galaxies  underlying  the  deduction  of  the 
two  tests  are  the  same  as  in  the  original  paper  of  1954;  they  are  reformulated  in 
the  following  section.  Admittedly,  these  assumptions  are  rather  specific  and  not 
very  likely  to  be  satisfied  exactly.  Later,  we  indicate  that  whenever  the  number 
of  observations  is  substantial  so  that  probability  limit  theorems  become  appli¬ 
cable,  the  domain  of  validity  of  the  tests  is  markedly  extended. 

We  remark  that  applications  of  the  virial  theorem  test  also  involve  assump¬ 
tions  concerning  the  systems  of  galaxies — sometimes  specific  (for  example,  [11]) 
and  sometimes  understood.  In  addition  to  the  assumptions  regarding  the  systems 
of  galaxies,  the  two  kinematic  tests  of  stability  depend  on  assumptions  of  a 
different  character,  concerned  with  the  collection  of  data.  Naturally,  assump¬ 
tions  of  this  kind  also  underly  the  virial  theorem  test  and,  indeed,  any  other 
test  based  on  radial  velocities  and  applied  to  individual  clusters  or  groups  of 
galaxies.  One  such  assumption  is  that  the  objects  treated  as  members  of  a  par¬ 
ticular  cluster  or  group  really  are  members  of  this  system,  rather  than  foreground 
or  background  galaxies.  Another  assumption  of  the  same  category  is  that  the 
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decision  as  to  whether  a  given  object  does  or  does  not  belong  to  the  particular 
system  of  galaxies  is  reached  without  reference  to  the  object’s  radial  velocity. 
These  two  important  assumptions  underlying  the  kinematic  tests  may  be  com¬ 
bined,  as  follows:  even  though  the  observations  relating  to  a  given  system  of 
galaxies  may  be  subjected  to  selection  for  brightness  or  for  position  or  for  some 
other  characteristic,  they  represent  a  fair  sample  of  the  system’s  membership  as 
far  as  radial  velocities  are  concerned. 

At  this  point  it  may  be  relevant  to  notice  that  while  the  validity  of  the  kine¬ 
matic  tests  is  not  harmed  in  any  way  by  selection  of  objects  for  their  brightness, 
this  probably  is  not  true  for  the  dynamical  tests.  The  reason  is  that  brightness 
generally  is  assumed  to  be  correlated  with  mass  and  thus  selection  for  brightness 
must  imply  selection  for  mass. 

Whether  the  data  relating  to  any  given  system  satisfy  the  combined  condition 
above  is  a  subjective  question.  In  order  to  give  the  reader  an  opportunity  to 
make  up  his  own  mind,  whenever  there  is  a  question  of  some  “outliers,”  which 
might  be  suspected  of  not  belonging  to  the  system  studied,  the  tests  for  stability 
were  applied  twice,  including  and  excluding  the  outliers.  In  fact,  the  summary 
test  was  applied  to  groups  of  galaxies  yet  again,  excluding  those  groups  that,  for 
one  reason  or  another,  appear  to  have  been  affected  by  selection  for  those  objects 
which  have  similar  radial  velocities. 

Quite  apart  from  the  above  precautions,  we  should  mention  a  particular  mech¬ 
anism  through  which  the  selection  of  galaxies  may  have  modified  the  distribution 
of  radial  velocities  and  invalidated  the  tests.  There  is  little  doubt  that,  particu¬ 
larly  for  distant  systems,  the  objects  selected  for  observation  are  those  with 
bright  apparent  magnitude,  including  the  dimming  due  to  redshift.  In  other 
words,  the  available  observational  material  does  not  include  galaxies  with  pecul¬ 
iar  velocities  away  from  the  observer  so  large  that  the  resulting  redshift  de¬ 
creases  their  apparent  brightness  to  the  point  of  making  them  difficult  to  observe. 

It  must  be  obvious  that  this  kind  of  selection  takes  place,  and  it  is  our  in¬ 
tention  to  study  it  in  another  paper.  For  the  present  we  note  that  the  values 
of  the  A-correction  to  the  apparent  magnitude  published  by  Humason,  Mayall 
and  Sandage  [13]  indicate  that  the  numerical  effect  of  the  selection  described  is 
likely  to  be  negligible. 

We  first  note  that  the  estimates  of  dispersion  of  radial  velocity  within  clusters 
are  of  the  general  order  of  1000  km/sec.  Thus  the  difference  between  the  radial 
velocities  of  two  cluster  members  selected  at  random  is  a  random  variable, 
probably  not  far  from  normal,  with  zero  mean  and  with  dispersion  of  about 
1400  km/sec.  This  difference  rarely  would  be  greater  than  3000  km/sec  and 
practically  never  greater  than  four  times  the  dispersion,  that  is,  5600  km/ sec. 
The  latter  figure  differs  but  little  from  the  radial  velocity  of  the  Perseus  Cluster 
(5433  km/sec)  for  which  Humason,  Mayall  and  Sandage  list  the  A-eorrection 
as  0.08  mag.  Even  allowing  for  the  nonlinearity  of  the  A-eorrection,  this  result 
indicates  that  the  differences  in  dimming  of  cluster  members  due  to  differences 
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in  their  radial  velocities  must  be  less  than  one-tenth  of  a  magnitude  and  can 
have  very  little  effect  on  the  decision  of  the  observer  to  include  a  given  object 
in  his  program. 

Whatever  the  case  may  be,  after  calling  the  reader’s  attention  to  the  theo¬ 
retical  possibility  of  bias,  in  the  remainder  of  the  present  paper  we  postulate 
that,  whatever  may  have  been  the  selection  of  galaxies  for  brightness,  this 
selection  does  not  affect  the  conditional  distribution  of  redshift.  What  we  assume, 
in  effect,  is  that  if  a  particular  galaxy  has  been  left  off  the  observational  program 
because  of  its  apparent  faintness,  then  this  galaxy  would  have  been  left  off  the 
program  also  if  it  were  brighter  by  one-tenth  of  one  magnitude. 

It  is  unfortunate  that  the  outcome  of  the  tests  depends  very  much  on  the 
objects  to  which  they  are  applied.  Thus,  for  example,  the  data  for  the  Coma 
Cluster  appear  to  include  four  “outliers,”  which,  for  various  reasons  explained 
in  the  sequel,  the  present  authors  suspect  do  not  really  belong  to  the  Coma 
Cluster.  However,  the  galaxies  in  question  are  located  in  the  general  area  of  the 
Coma  Cluster  and,  in  spite  of  indications  to  the  contrary,  may  conceivably 
belong  to  the  cluster.  If  these  outliers  are  accepted  as  members  of  the  cluster, 
the  test  leaves  little  doubt  that  the  cluster  is  expanding.  On  the  other  hand,  if 
the  outliers  are  omitted,  the  same  test  fails  to  show  expansion  and,  if  anything, 
indicates  that  the  Coma  Cluster  is  contracting.  From  the  point  of  view  of  the 
hypothesis  that  contraction  or  expansion  is  a  relatively  simple  phenomenon 
applying  to  the  cluster  as  a  whole,  these  contradictory  conclusions  certainly  are 
disconcerting.  However,  it  is  notable  that  qualitatively  they  agree  with  the 
following  conclusion  of  van  Albada  [14]  suggested  earlier  by  Oort  [10] :  “At  a 
later  stage,  relative  contraction  of  the  central  regions  (of  a  cluster)  becomes 
absolute  contraction,  while  the  intermediate  and  outer  regions  are  still  expand¬ 
ing.”  There  is,  then,  the  possibility  that  the  suspected  “outliers”  in  the  Coma 
Cluster  actually  are  members  of  this  system  but  are  located  in  the  “outer 
region”  showing  a  rapid  expansion.  At  the  same  time,  the  central  region  of  the 
cluster  may  be  contracting. 

A  somewhat  similar  dependence  of  the  conclusions  on  the  material  included 
in  the  calculations  prevails  in  the  study  of  groups  of  galaxies.  The  two  summary 
tests  applied  to  all  groups  for  which  we  have  data  fail  to  show  evidence  of  in¬ 
stability.  On  the  other  hand,  if  those  groups  for  which  there  is  a  suspicion  of 
selection  for  similarity  in  redshift  are  omitted,  then  the  same  tests  yield  a  strong 
indication  of  expansion. 

The  above  are  examples  of  objective  findings  presented  below.  Our  subjective 
opinion  is  that  the  hypothesis  that  at  least  some  systems  of  galaxies  are  unstable 
cannot  be  dismissed  lightly  and  deserves  further  study. 

Certain  remarks  by  Dr.  R.  Minkowski  made  at  the  Symposium  during  the 
discussion  of  the  present  paper,  combined  with  some  further  analysis  presented 
below,  lead  us  to  another  conclusion.  This  is  that,  if  the  results  of  the  kinematic 
tests  indicating  instability  really  reflect  instability  of  systems  of  galaxies,  rather 
than  an  arbitrary  selection  of  data,  then  the  same  results  are  evidence  of  the 
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presence  of  absorbing  matter  within  the  systems  concerned.  The  study  of  this 
question  is  in  progress. 

Throughout  this  paper,  the  term  “system  of  galaxies”  is  used  as  a  synonym 
of  the  term  “cluster”  appearing  in  the  theory  of  the  spatial  distribution  of 
galaxies  [1],  [15].  This  change  in  terminology  is  adopted  in  deference  to  the 
distinction  made  in  some  sections  of  astronomical  literature  in  which  the  word 
cluster  is  used  to  denote  a  system  of  a  large  number  of  galaxies,  while  a  system 
of  relatively  few  members  is  described  as  a  group. 

Although  the  methods  used  in  this  paper  are  applied  solely  to  systems  of 
galaxies,  the  theoretical  results  described  in  the  next  sections  could  as  well  be 
adapted  to  stellar  systems. 


2.  Kinematic  test  of  the  stability  of  several  systems  of  galaxies 

against  the  alternative  that  all  these  systems  expand 

or  contract  at  the  same  rate 

In  this  section  we  recapitulate  and  slightly  reformulate  the  assumptions  under¬ 
lying  the  deduction  of  the  locally  most  powerful  test  for  stability  of  a  system  as 
given  in  [1]  and  indicate  an  easy  generalization  of  this  test  so  as  to  refer  to  any 
number  s  of  systems  assumed  to  have  certain  similar  properties. 

The  basic  hypotheses  underlying  the  test  of  stability  are  those  adopted  in  the 
Neyman-Scott  theory  of  the  spatial  distribution  of  galaxies.  Consider  an  orthog¬ 
onal  system  of  coordinates  with  its  origin  at  the  observer.  Consider  a  system 
with  its  center  at  u  =  {uh  u2,  u3)  and  let  D  be  the  distance  of  this  center  from 
the  observer.  Our  hypotheses  are : 

(i)  Given  u,  the  coordinates  X  =  (Xi,  X2,  X3)  of  any  member  of  the  system 
are  random  variables  with  their  conditional  probability  density  /(r?)  depending 
only  on  the  distance  t?  from  the  center  of  the  system. 

In  the  applications  of  the  theory,  whenever  we  need  to  deal  with  specific 
functions,  it  is  assumed  that 

(1)  fM  ~ 

(ii)  Given  the  positions  uh  u2,  •  •  •  ,  us  of  the  centers  of  s  systems,  the  positions 
of  the  objects  belonging  to  these  systems  are  mutually  independent  and  follow 
the  same  distribution. 

(iii)  The  absolute  magnitude  M  of  a  cluster  member  is  a  random  variable 
independent  of  all  the  other  random  variables  considered  and  following  a  normal 
distribution  with  mean  M0  and  variance  oh- 

(iv)  The  present  study  was  conducted  from  the  point  of  view  of  the  hypoth¬ 
esis  that  the  amount  of  absorbing  matter  in  space,  whether  within  the  systems 
considered  or  outside,  is  negligible.  However,  in  order  to  explain  certain  details 
of  the  numerical  findings,  we  feel  in  need  of  the  hypothesis  that  the  systems 


192  FOURTH  BERKELEY  SYMPOSIUM:  LOVASICH,  MAYALL,  NEYMAN,  SCOTT 

can  include  some  dark  matter.  The  actual  distribution  of  this  matter  is  likely  to 
be  complex.  However,  it  is  hoped  that  the  following  simple  hypothesis  may  serve 
as  a  first  approximation. 

Consider  a  system  with  its  center  at  distance  D  from  the  observer.  Consider  a 
particular  member  of  this  system  and  denote  by  £  its  distance  from  the  observer 
and  by  <p  its  angular  separation  from  the  center  of  the  system.  Let  C  denote  the 
dimming  of  an  object  located  exactly  at  the  center  of  the  cluster,  that  is,  when 
£  =  D  and  <p  —  0.  We  assume  that,  for  small  <p,  the  same  number  C  will  represent 
the  dimming  of  all  the  cluster  members  for  which  £  =  D  cos  <p.  Generally,  for 
small  <p  and  arbitrary  £,  we  assume  that  the  apparent  magnitude  m  (corrected 
for  redshift)  is  connected  with  the  absolute  magnitude  of  the  object  by  the 
equation 

(2)  m  —  M  —  5  +  5  logi0£  +  C  +  r(£  —  D  cos  <p), 

where  r  is  a  constant  to  be  described  as  the  absorption  coefficient.  The  value  of 
C  does  not  play  any  role  in  our  calculations.  As  mentioned,  the  principal  com¬ 
putations  were  performed  assuming  r  =  0.  However,  certain  details  indicate 
that  t  may  be  a  positive  number  and  some  of  the  computations  were  repeated. 

(v)  We  now  consider  the  velocities  of  cluster  members  with  respect  to  the 
center  of  the  system  to  which  they  belong. 

The  velocity  vector  £2  of  a  cluster  member  located  at  X  is  the  resultant  of 
two  vectors  £2i  and  fl2,  the  first  of  which  is  random  and  the  second  nonrandom. 
Given  the  positions  of  several  cluster  members,  the  random  components  are 
mutually  independent  and  identically  distributed.  More  specifically,  the  distri¬ 
bution  of  the  component  12i  =  (Xh  X2,  X3)  is  spherical  normal,  with  variance 
<7?  and  with  E(Xj)  =  uj,  where  j  =  1,  2,  3  and  u  =  (uh  u2,  u3)  is  the  recession 
velocity  of  the  center  of  the  system  to  which  the  given  object  belongs.  As  to 
the  nonrandom  vector  tt2,  it  represents  the  velocity  of  the  Hubble-type  expan¬ 
sion  (or  contraction)  of  the  system.  This  vector  is  assumed  to  be  directed  away 
from  the  center  of  the  system  and  its  magnitude  is  assumed  to  be  proportional 
to  the  distance  between  the  center  and  the  cluster  member  in  question.  Thus, 
for  a  member  with  fixed  coordinates  X\,  x2,  x3,  the  components  of  fi2  are 

(3)  6(xi  -  Mi),  d(x2  -  m2),  0(x3  -  m3), 

w’here  6  is  the  factor  of  proportionality,  analogous  to  the  Hubble  constant.  If 
0  =  0  then  the  system  is  stable,  while  9  >  0  means  expansion  and  9  <  0  means 
contraction.  Figure  1  illustrates  the  situation. 

(vi)  It  is  assumed  that  the  radial  velocity  of  the  center  of  each  system  is 
exactly  proportional  to  the  distance  of  this  center  from  the  observer,  so  that 
the  components  of  this  velocity, 

(4)  iij  =  Hujf  j  =  1,  2,  3, 

where  H  is  the  Hubble  constant. 

The  variables  that  are  observable  for  each  object  are  its  radial  velocity  p,  its 
apparent  magnitude  m,  which  we  shall  assume  to  have  been  corrected  for  red- 
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shift,  and  its  angular  distance  <p  from  the  center  of  the  system  considered.  Of 
these  three  random  variables  we  shall  assume  that  the  data  collected  for  any 
system  are  subject  to  selection  for  m  and,  possibly,  for  <p.  On  the  other  hand, 
we  assume  that  there  is  no  selection  for  p.  The  test  of  stability  of  the  system  is, 
essentially,  a  test  of  the  hypothesis  that  the  observed  values  of  p  conform  to  the 


Random  and  nonrandom  velocity  components. 

For  illustrative  purposes,  two  cluster  members  are  assumed  to 
be  at  the  same  distance  i\'  from  the  center  of  the  system. 

They  are  denoted  G'.  Two  other  cluster  members,  denoted  G", 
are  assumed  to  have  distance  from  the  center  rj"  =  2?j'.  Thus 
the  magnitudes  of  the  corresponding  vectors  Q>  are  in  the  same  relation. 

conditional  distribution  of  this  variable  implied  by  the  assumptions  (i)  to  (vi) 
in  combination  with  the  hypothesis  of  no  expansion  0  =  0. 

Geometrical  considerations,  combined  with  the  hypotheses  regarding  the  ve¬ 
locity  vector  ft,  imply 

(5)  p  =  \  [  £  xtx,  +  e(p  -  £ 

f  L  i=i  i=i  J 

=  v  cos  +  HD  cos  <p  +  0(£  —  D  cos  <p), 

where  v  =  |fti|  is  the  random  component  of  the  peculiar  velocity  and  is  the 
angle  between  fti  and  the  line  of  sight.  Formula  (5)  implies  that,  for  fixed  X 
and  u,  the  conditional  distribution  of  p  is  normal  with  variance  af  and  with 
expectation 

(6)  E{p\  X,  u)  =  HD  cos  <p  +  0(£  —  D  cos  <p). 

Consider  s  systems  and  assume  that  the  distribution  of  <M,  of  Xi  —  Ui,  and 
of  p  are  identical  from  one  system  to  the  next.  Also  we  assume  in  this  section 
that  the  rate  of  expansion  8  is  the  same  for  each  system.  The  subscripts  i  or  ij 
attached  to  various  symbols  will  indicate  the  fth  system  and  the  jth  member  in 
the  tth  system,  respectively.  We  shall  have  i  =  1,  2,  *  •  •  ,  s  and,  for  a  fixed  i, 
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the  subscript  j  =  1,  2,  •  •  •  ,  n,-,  where  n,  stands  for  the  number  of  objects  in 
the  ith  system  for  which  the  necessary  data  are  available.  Our  purpose  is  to 
deduce  the  joint  conditional  probability  density,  say  p,  of  all  the  ptj,  given  the 
values  of  all  the  and  Because  of  the  postulated  independence,  p  is  the 
product  of  the  conditional  densities  of  the  pa  given  and  ho¬ 
using  familiar  notation  and  suppressing  the  subscripts  i,  j  for  the  time  being, 
we  have 


(7) 


Pp\m,<p 


Pm,<p,p 

Pm,<p 


In  order  to  simplify  the  formulas,  rather  than  work  with  the  variables  m,  p  and 
£,  it  will  be  convenient  to  deal  with  an  equivalent  system  defined  by 


ju  =  m  -f-  5  —  5  logioCD  cos  <p)  —  C  —  M0, 


(8) 


r  =  p  —  HD  cos  <p, 

_  £  —  D  cos  <p 
a 


The  adopted  hypotheses  then  imply  that 

p-D* sinV/2<rs 

(9)  Pp,t  =  - D 2  cos2  ip  sin  <p  e~J2/2(l  +  qt )2, 

a2V2x 

where  q  =  a/D  cos  <p  and  the  region  of  variation  of  t  extends  from  —q1  to 
infinity.  Given  t  and  <p,  the  conditional  density  of  p  is  then 

(10)  p^.t  =  - K=  exp  {  -  7TT  [m  -  a  log(l  +  qt)  -  rat]2  j, 

o-,wv2x  l  Z(Jm  J 

where  the  logarithm  is  taken  to  the  base  e  and  a  =  5  logi0e  =  2.1715.  Finally, 
the  conditional  probability  density  of  r,  given  tp  and  t,  is 


(11)  VrW.t  =  e-fr-“>V2. 

0"1  V  Z7T 


Familiar  formulas  give 


(12) 

Pn,<p  j—g-i  PvlvA  dt, 

(13) 

Pn,<p,r  J  V<p,t  Pfi\<p,t  Pr\<p,t  dt, 

and  finally 

J ( r ,  m,  <p) 

(14) 

p*w  Hf,  v)  ’ 

with 

(15)  J  (r,  m,  <p) 

=  /_”.,(!  +  Qt)2  exp  | 

1  I"  p  +  I'm  -  a  log(l  +  qt)  ~ 

2  |_  0Af 

(r  —  6  at)2  T) 

<A  Jj 


dt, 
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(16)  I(n,  <?) 

=  ^  (1  +  qty  exp  {  -\ [  f2  +  dL 

The  joint  probability  density  p  is  then 

(17)  P  =  II  n  Pm^.Pir 

1  =  1 j= 1 

Formulas  (14)  to  (16)  indicate  that,  in  general,  the  conditional  distribution 
of  r  given  p  and  <p  is  not  normal;  this  creates  the  problem  of  deducing  a  test  of 
the  hypothesis  0  =  0  satisfying  some  optimum  condition.  However,  it  will  be 
noted  that,  if  the  hypothesis  tested  is  true,  then,  for  fixed  p  and  <p,  the  variable  r 
is  normally  distributed  with  mean  zero  and  variance  oi.  This  is  an  important 
detail  which  simplifies  the  problem  considerably.  The  test  to  be  deduced  has  the 
property  of  being  locally  most  powerful  unbiased. 

Repeating  the  reasoning  explained  in  [1],  we  find  that  the  locally  best  test  of 
the  hypothesis  that  0  =  0  reduces  to  the  rule  of  rejecting  this  hypothesis  when¬ 
ever  the  absolute  value  of  the  criterion 

(18)  T  =  ~ 

<T$ 

exceeds  the  tabled  value  of  Student’s  t  corresponding  to  the  chosen  level  of 
significance  a  and  to  the  number  of  degrees  of  freedom  £wt  —  s  —  1.  Here  0 
stands  for  the  least  squares  estimate  of  0, 

X  XI  (-d- ij  -T 1')  Pi]  COS  (fij 

(19)  9  =  E  E  (Ai,  -  A,.)2cosVi,- 

1  j 

with 


(20) 


\  _  G  ^  \tjj *Pij)  _  &  \S  jj 


COS  <Pij 


COS  (pij 


As  to  Ai.  and  a  similar  symbol  pi.  which  we  shall  need  later,  they  are  given  by 

A.  ij  COS  tp ij  ^  v  Pij  COS  (p ij 

(21)  A,.  =  4= - - - -  =  -(= - ? - • 

2.  cos2  va  2-  cos2  tpij 

}  3 

The  quantities  A  a  can  be  computed  by  numerical  integration  using  (9)  and  (10). 
These  formulas  depend  on  the  distance  D,  to  the  cluster.  In  computations  these 
distances  must  be  replaced  by  their  estimates  pi-H~l.  The  numerical  integration 
involved  is  tedious  and  further  below  we  indicate  convenient  approximations 
based  on  the  fact  that,  ordinarily,  the  quotients  q,j  =  a/D,  cos  <pa  and  the  angles 
(pij  are  small. 

The  quantity  &§  in  the  denominator  of  (18)  is  the  estimate  of  the  standard 
error  of  0.  The  variance  of  0  is  given  by 
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(22) 


L  L  (Ai)  -  Ai.)2  cos2  <pij 

i  } 


In  order  to  obtain  its  estimate  to  be  used  in  (18),  the  value  of  a\  in  (22)  is  re¬ 
placed  by 


(23) 


2 Z  \Pij 

* 2  /  i _ 

o  1  — 


X  Ui  —  s  —  1 


which  is  the  pooled  estimate  of  the  conditional  variance  <j\  of  p. 

In  the  above  formulas,  the  summations  for  i  always  extend  over  all  the  systems 
of  galaxies  included  in  the  test,  that  is,  over  i  =  1,  2,  •  •  •  ,  s.  The  summations  for 
j  extend  over  j  =  1,  2,  •  •  •  ,  n,-.  When  we  put  s  =  1  the  formulas  of  this  section 
become  equivalent  to  those  in  the  original  publication  [1].  Whatever  the  value 
of  s,  when  the  hypothesis  tested  is  true,  the  T-test  will  reject  it  falsely  with 
probability  exactly  equal  to  a,  irrespective  of  the  values  of  the  D%  and  of  a i, 
and  no  matter  how  the  objects  were  selected  for  magnitude  and  position. 

It  may  be  useful  to  point  out  the  relation  between  the  T-test  just  described 
and  the  classical  t- test.  Operationally,  the  two  coincide.  However,  the  /-test  has 
the  property  of  being  the  uniformly  most  powerful  test  corresponding  to  the 
assumption  that,  whether  the  hypothesis  tested  is  true  or  not,  the  observable 
random  variables  are  normally  distributed  with  a  fixed  variance.  Because  of  the 
nonnormality  of  the  distribution  of  r(and  therefore  of  p),  the  T-test  is  not  uni¬ 
formly  most  powerful  (no  such  test  exists)  and  it  has  a  power  function  different 
from  that  of  t.  However,  when  the  number  of  observations  is  large  and  0  is  small, 
the  power  function  of  T  is  approximately  equal  to  that  of  t. 


3.  Kinematic  test  of  the  stability  of  several  systems  of  galaxies 
against  the  alternative  that  they  may  be  expanding 
or  contracting  each  at  a  possibly  different  rate 

In  this  section  we  contemplate  the  possibility  that  among  the  systems  studied 
there  may  be  some  that  are  in  fact  stable,  while  some  others  may  be  expanding 
each  at  its  own  rate,  and  still  others  may  be  contracting.  In  other  words,  con¬ 
trary  to  what  was  assumed  in  the  previous  section,  we  now  assume  that  to  the 
zth  system  of  objects  there  corresponds  its  own  rate  of  expansion  0t  which  need 
not  be  equal  to  the  rate  Qj  of  the  j'th  system.  The  hypothesis  to  test  is  now 

(24)  0i  =  02  =  •  •  •  =  0,  =  0. 

All  other  assumptions  used  in  section  2  are  maintained.  The  analogy  of  the 
T-test  to  the  classical  t- test  suggests  the  possibility  of  extending  the  T-test  in 
the  same  manner  in  which  the  classical  Z- test  is  extended  for  testing  the  general 
linear  hypothesis.  As  is  well  known,  this  involves  the  use  of  the  familiar  T-dis- 
tribution. 
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Considering  only  one  system  of  objects,  say  the  tth,  and  using  formula  (19)  we 
obtain  the  least  squares  estimate  of  d„ 


(25) 


]X  (d-  ij  Aj.)pjj  COS  <p ; 

9l  E  (A,j  -  Ai.y  cos2  (pi3 


The  variance  of  0*  is  obtained  from  a  similar  particularization  of  (22),  namely, 


(26) 


<r“ 

0* 


X  (Aij  —  A,-.)2  cos2  <pij 


However,  because  of  the  admitted  possibility  that  each  system  of  galaxies  may 
be  unstable  at  a  different  rate  0,,  the  pooled  estimate  of  a'{  is  no  longer  given  by 
(23)  but  by  the  formula 


(27) 


aj  = 


E  E  (pij  -  pi  )-  -  E  0i  E  (*4 o 

-i _ i _ * _ i _ 

X  m  -  2  s 


Ai)pij  cos  <p,j 


with  the  number  of  degrees  of  freedom  Ew»  —  2 s. 

If  the  hypothesis  tested  is  true,  then  each  product 

X  (A  ij  A{.)pij  COS  <Pij 

(28)  h  IjE  (A«  -  Ai.y  cos2  TOJ»  =  j-£  (AiS  -Ai.y  cos2^J'2 

is  a  normal  variable  with  expectation  zero  and  variance  a\.  It  follows  that  in 
this  case  the  variable,  say, 

E  0?  E  (A ij  -  Ai.y  cos2  <prj 

J. _ i _ _ _ 

S  ai 

follows  exactly  the  l '-distribution  with  degrees  of  freedom  fx  =  s  and 
$2  —  n  i  —  2s.  Thus,  the  test  of  the  hypothesis  0i  =  02  =  •  •  •  =  0.s  =  0  reduces 
to  the  rule  of  rejecting  this  hypothesis  whenever  the  computed  value  of  Z  exceeds 
the  tabled  value  of  F  corresponding  to  the  indicated  degrees  of  freedom  and  to 
the  chosen  level  of  significance  a.  When  the  hypothesis  test  is  true,  it  will  be 
rejected  falsely  with  probability  exactly  equal  to  a,  irrespective  of  the  values  of 
the  Di  and  of  ah  and  no  matter  how'  the  objects  were  selected  for  magnitude 
and  position. 


4.  Approximate  formula  for  the  symbol  A 

According  to  definition  (20)  and  using  (9)  and  (10),  we  have 

[ t(  1  +  qty  exp(  {£2  +  0  —  a  log(l  +  qt)  - 

(30)  A  cos  ^  —  a  - - — - 

/  (1  +  qt)2  exp(  —  \  (i f 2  +  [m  —  a  log(l  +  qt)  - 


Toty I a\f))  dt 
rat]2 / aiir)  )  dt 
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where  both  integrations  extend  from  —  q~x  to  -f  and  the  subscripts  ij  have 
been  suppressed  for  brevity.  Using  the  value  of  a  in  the  range  105  to  10®  as 
estimated  in  [16]  and  the  estimates  of  the  distances  Di  of  the  various  systems 
of  galaxies  involved  in  the  empirical  study  presented  below,  we  find  that,  except 
for  the  closest  group,  q a  =  a/Dl  cos  <pa  is  always  smaller  than  0.05.  This  fact 
indicates  that  the  integrals  in  (30)  depend  mainly  on  the  properties  of  the  inte¬ 
grands  in  the  interval  |l|  <4.  Thus,  with  the  danger  of  committing  only  a  small 
error,  we  may  consider  these  integrals  as  extending  from  —  oo  to  +<*>  and  also 
we  may  use  the  approximation 


(31)  log(l  +  qt)  =  qt. 

Adopting  these  approximations  and  performing  easy  transformations,  we  write 


(32) 


A  cos  <p  =.  a 


J 

f  £(1  -f-  qt)2  exp 

Cm  +  XV2  .  v. 

-  02  \t~  W 

Zcm 

dt 

a  ~ 

J 

f  (1  +  Qt)2  exp 

cm  +  XV  ~ 

.  -  2,5,  ('  “  (o)  J 

dt 

where,  for  the  sake  of  brevity, 

(33) 


X  =  ^+r, 
a 


(34) 


Xo~M 

om  +  XV 


and  the  two  integrals  extend  from  —  «>  to  +°o.  Simple  calculations  now  give 


(35) 


A  cos  <p 


a  to 


+ 


2(7(1  +  qto)alf 


(1  +  qto)2(<TM  +  XV)  +  q2CMj 

Already,  this  formula  for  A  cos  is  easy  to  compute.  However,  if  one  decides 
to  neglect  the  terms  in  q 2  and  also  to  substitute  unity  for  cos  <p,  which  introduces 
errors  of  the  order  of  a  few  percent  at  most,  then 

_  q-(Xo-m  +  2qalr) 

<r\i  +  XV 


(36) 


A 


r  r?  1 

,  .  0  Cm  __ 

<’P„  +  2vm  +  xvJ  - 


Numerical  computations,  involving  the  evaluation  of  integrals  by  quadrature 
formulas,  indicated  that,  for  all  but  the  closest  systems,  formula  (36)  provides 
a  very  good  approximation.  For  the  closest  systems  the  approximation  to  indi¬ 
vidual  values  of  the  A  a  is  somewhat  less  satisfactory.  However,  since  the  A  i} 
appear  in  all  statistics  as  weights,  the  final  effect  of  the  inaccuracies  is  negligible. 
For  example,  in  the  Ursa  Major  Cloud,  where  <ptj  extends  to  15°,  the  error  in 
the  estimate  of  di  is  only  five  percent. 

Returning  now  to  the  subscripts  ij,  formula  (36)  yields 


(37)  A  u  -  A  i.  =  2  ( ml]  -  mi.), 

cm  “r  ^iC 

with 

(38)  X;  =  r  +  JJ- 
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Before  proceeding  further,  we  recall  that  the  symbol  m  is  used  here  to  denote 
the  apparent  magnitude  of  a  galaxy  corrected  for  redshift.  On  the  other  hand,  the 
apparent  magnitudes  usually  are  published  without  this  correction,  which  is 
often  called  the  K-correction.  The  formulas  given  above  could  be  used  directly 
after  each  of  the  apparent  magnitudes  had  been  corrected.  However,  a  con¬ 
siderable  amount  of  labor  can  be  saved  by  noticing  that  over  a  very  substantial 
range  of  redshift  the  ^-corrections  are  very  small  and  are  approximately  pro¬ 
portional  to  p.  In  fact,  using  the  ^-corrections  published  by  Humason,  Mayall 
and  Sandage  [13]  for  twelve  clusters  ranging  in  radial  velocities  up  to  23,365 
km/sec,  it  was  found  that,  with  only  trivial  errors,  we  can  put 

(39)  K(p )  =  xp  with  x  =  1-55  X  10~5. 

Now  if,  contrary  to  the  notation  adopted  earlier,  m  denotes  the  apparent  mag¬ 
nitude  not  corrected  for  redshift,  then  expression  (37)  for  A  a  —  A,,  must  be 
replaced  by 

X  <72 

(40)  Aij  A{.  =  ~ 2  i  \2  9.  [(^L'i  m,.)  x(po'  P**)]* 

<?m  -r  Ai<r 

Now  denote  by  $?(p)  the  sample  variance  of  p,  by  Sfm,  p)  the  sample  covariance 
of  p  and  the  uncorrected  apparent  magnitude  and,  finally,  by  Sf(w)  the  sample 
variance  of  the  uncorrected  apparent  magnitude,  all  computed  from  the  data  on 
ni  objects  within  the  ith  system,  so  that 

£i(p)  =  ]_  1  L  (po-  -  p.-)2, 

(41)  Si(m,  p)  =  w  l_  ■  Y  (mu  -  mi.)(Pij  -  Pi), 

fli  1  j 

S i(m)  =  — “j-  Y  (ma  ~  m,.)2. 
n%  i  j 

For  purposes  of  computing,  the  main  elements  in  the  formulas  connected  with 
the  two  tests  for  stability  are  easy  to  express  in  terms  of  (41).  We  have 


(42) 

with 

(43) 


X.  (An  -  a,.)*  cos2  ip,j  N  [A^-J  (m  -  l)ci(m), 

X  a*2 

Y*  (-4. i>  Ai-^pij  cos  <pij  •  2  i  \2  2  (a*  1  )<x,(?w,  p), 

j  ~r  Ai<r 

<f\(m)  =  S\(m)  -  2 xSi(m,  p)  +  x2S2t(p), 

<n(m,  p)  =  Si(m,  p)  -  xSt(p) 


representing  the  sample  variance  of  the  apparent  magnitude  corrected  for  red¬ 
shift  and  its  sample  covariance  with  redshift,  respectively.  Now  the  expressions 
in  (42)  can  be  substituted  into  the  formulas  giving  the  estimates  6,  0»,  &\,  <r|  and 
so  forth  used  in  the  T-test  and  in  the  Z-test.  Using  (19)  and  (42)  we  have 
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(44) 

Also, 

(45) 


S  ( —  l)X,-(Ti(m,  p)/ (<&  +  XV) 

0  *  O*  E  (»<  -  1)  XW(m)/(A  + 

1 

0.2 

N  <r4  £  (n<  -  l)XW(m)/(A  +  *VF 

i 


and  so  forth. 

When  considering  the  formulas  above,  it  is  important  to  be  clear  how  far  their 
values  depend  on  the  data  involved  in  the  study  and,  on  the  other  hand,  how 
far  they  are  influenced  by  constants  taken  from  other  sources  and  possibly  sub¬ 
ject  to  doubt.  Such  constants  are  x,  &m,  <t,  t  and  the  Hubble  constant  H. 

The  values  of  the  ^-correction  are  deduced  using  certain  hypotheses  and 
approximations,  and  depend  on  the  type  of  galaxies  involved.  However,  unless 
the  theory  behind  them  is  entirely  wrong,  the  order  of  magnitude  of  the  figures 
given  by  Humason,  Mayall  and  Sandage  may  be  expected  to  be  correct.  The 
values  of  the  expressions  (44),  (45)  and  others  will  not  change  much  if  x  =  1-55 
X  10-6  were  to  be  multiplied  or  divided  by  a  substantial  factor.  Thus,  we  do  not 
anticipate  any  difficulty  from  the  value  of  x- 

The  other  constants  <r,  cM,  r  and  H  are  much  more  troublesome.  We  can  write 
H  =  10 ~*h  km/sec/pc  where  h  is  a  factor  between  0.5  and  1.5,  depending  on  the 
distance  scale  one  wants  to  adopt.  Then,  using  the  assumption  in  (4),  we  estimate 
the  distance  to  the  ith  cluster  by  Z>»  =  H~xpi.  =  104p,-./X  parsecs.  The  value  of 
the  dispersion  a  also  depends  on  H,  namely,  a  =  9/h  where  ?  is  of  the  order  of 
7  X  106  parsecs.  On  the  other  hand,  the  dispersion  <rM  does  not  depend  on  H  ;  the 
value  <tm  =1.25  mag  used  below  represents  a  compromise  between  several  cur¬ 
rent  estimates.  As  to  r,  there  is  an  estimate  by  Zwicky  [9]  to  the  effect  that 
galaxies  located  behind  a  large  cluster  have  their  apparent  magnitudes  increased 
by  about  0.7  mag  due  to  absorbing  matter  within  the  cluster.  Assuming  that 
the  radius  of  a  cluster  corresponds  to  about  3.5 <r,  Zwicky ’s  estimate  yields 
<tt  —  0.1  mag. 

Taking  all  these  values  into  account,  we  notice  that 


is  independent  of  h  and,  for  not  too  small  p,.,  relatively  unimportant  compared 
to  &m .  It  follows  that  even  substantial  errors  in  r  will  not  affect  the  denominators 
om  +  X?<r2  appearing  in  (44)  and  (45).  On  the  other  hand,  the  same  approximate 
calculations  indicate  that  if  the  adopted  value  of  is  an  underestimate  or  an 
overestimate  of  the  true  value  by  a  factor  then  the  estimate  0  is  multiplied  by 
the  same  factor.  Finally  the  factor  <r2  appearing  in  the  denominator  of  (44) 
indicates  that  the  estimate  0  is  roughly  proportional  to  h2. 

To  summarize,  we  may  say  that  the  estimates  0  and  0<  are  nearly  proportional 
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to  a\t/ a1  and  that  is  nearly  proportional  to  the  square  of  the  same  ratio.  This 
circumstance  must  be  taken  into  account  when  examining  the  numerical  esti¬ 
mates  obtained. 

Should  it  happen  that  the  outcomes  of  the  T  and  the  Z-tests  depend  on  the 
assumed  values  of  the  above  parameters  to  any  comparable  extent,  this  would 
constitute  a  considerable  weakness  in  the  tests.  Luckily,  this  is  not  the  case. 
The  two  tests  depend  essentially  on  quotients  of  the  type  9/&g  and  it  so  happens 
that  the  factors  c^m/g1  occur  in  both  the  numerator  and  the  denominator  and 
cancel.  It  follows  then  that,  while  the  estimates  of  the  expansion  factors  do 
depend  on  the  uncertain  quantities  <tm  and  a,  the  outcomes  of  the  kinematic 
tests  are  practically  independent  of  these  parameters  and  depend  only  on  the 
data.  As  far  as  the  validity  of  the  tests  is  concerned,  the  greatest  danger  comes 
from  the  possible  selection  on  p. 

The  fact  that  the  outcomes  of  the  kinematic  tests  do  not  depend  on  the  values 
of  the  assumed  parameters  is  not  surprising  when  we  consider  the  intuitive 
interpretation  of  the  tests.  Essentially,  they  measure  the  correlation  between 
Aij  and  p,y  (using  the  approximation  for  An,  the  correlation  between  ra.y  and  p„) 
and  thus  are  independent  of  the  value  of  any  scale  or  location  parameter.  The 
heuristic  explanation  of  Ti}  for  example,  is  that  if  there  is  expansion  of  the  ith 
cluster  then  the  objects  on  the  front  of  the  cluster  (those  with  negative  Ai}) 
will  have  a  component  S22  towards  the  observer  and  thus  will  tend  to  have  smaller 
Pn  than  average,  whereas  those  objects  on  the  back  side  of  the  cluster  will  have 
positive  A  ,y  and  will  tend  to  have  larger  p,y  than  the  average.  Hence,  expansion 
will  tend  to  produce  large  positive  values  of  TV 

One  more  remark  is  necessary.  An  inspection  of  the  formulas  (14)  to  (16) 
indicates  that  But  is  the  conditional  expectation  of  r  =  p  —  D  cos  <p  whether  the 
hypothesis  tested  is  true  or  not.  This  circumstance  is  the  reason  why  the  esti¬ 
mate  6  is  a  good  test  criterion :  if,  relative  to  its  standard  error,  6  is  too  far  from 
zero  then  this  is  a  good  indication  that  the  true  9  differs  from  zero.  On  the  other 
hand,  it  is  necessary  to  realize  that  the  property  that  <rf  is  an  unbiased  estimate 
of  <7i,  the  conditional  variance  of  p  given  the  position  of  the  galaxy  within  the 
cluster,  has  a  “local”  character  and  is  subject  to  the  assumption  that  0  =  0. 
This  is  satisfactory  unless  one  adopts  the  attitude  that  0=^0.  When  MO  then 
<ti  becomes  the  estimate  of  the  conditional  variance  of  p  given  <p  and  m  and 
includes  the  variability  of  p  due  to  the  variability  in  the  distance  £  with  which 
p  is  correlated  and  which,  for  any  fixed  m,  must  be  considerable.  For  this  reason, 
the  formulas  for  the  standard  errors  of  the  estimates  9  are  not  appropriate  when 
it  is  admitted  that  9^0.  The  problem  of  estimating  the  standard  errors  in  this 
case  must  be  relegated  to  another  study. 

The  practical  application  of  the  two  tests,  in  their  approximate  form  indicated 
in  (44)  and  (45)  involve  the  following  preliminary  steps: 

(i)  Using  the  magnitudes,  not  corrected  for  redshift,  and  the  radial  velocities, 
compute  for  each  system  the  sample  variances  S^im)  and  $<(p)  and  the  covariance 
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Si(m,  p )  as  indicated  in  (41).  These  statistics  are  interesting  for  their  own  sake. 

(ii)  Using  the  estimates  Z),  =  H~xpi .  and  an  assumed  value  of  r,  compute 
X,  =  r  +  a/D{. 

(iii)  Using  an  assumed  value  of  x,  compute  < Aim)  and  <Xi(m,  p),  the  corrected 
variance  and  covariance,  from  (43). 

(iv)  Assuming  values  for  aM  and  for  a,  and  using  the  approximations  in  (42), 
proceed  to  calculate  7\,  T  and/or  Z. 

5.  Extended  validity  of  the  kinematic  tests 

As  noted  at  the  outset,  the  hypotheses  (i)  to  (vi),  used  to  deduce  tests  of 
stability  of  systems  of  galaxies,  are  not  likely  to  be  satisfied  exactly  by  actual 
systems  of  galaxies.  For  example,  to  mention  only  a  few  points,  it  is  likely  that 
the  distribution  of  positions  of  the  particular  galaxies  within  a  system  does  not 
correspond  to  formula  (1).  The  same  applies  to  the  distribution  of  the  random 
component  Hi  of  the  velocity  vector  of  a  galaxy.  Also  the  postulated  independence 
of  the  absolute  magnitude  M  and  of  the  random  velocity  component  Hi  is  subject 
to  doubt.  The  more  massive  galaxies  are  likely  to  be  bright  and  to  have  some¬ 
what  smaller  peculiar  velocities.  Furthermore,  M  need  not  be  normally  dis¬ 
tributed.  Finally,  the  vector  Hi  probably  depends  on  the  position  of  the  galaxy 
within  the  system. 

The  purpose  of  this  section  is  to  indicate  that,  while  the  optimality  of  the 
tests  deduced  does  depend  on  the  hypotheses  enumerated  above,  the  two  tests 
remain  approximately  valid  in  much  broader  conditions.  We  combine  these 
conditions  in  the  following  definition  of  a  “stable  system.” 

With  reference  to  figure  2,  let  C  be  the  center  of  a  system  of  galaxies  and  P  an 
arbitrary  point  in  space.  Draw  an  arbitrary  line  L  through  P  and  select  a  posi- 


Figure  2 

Illustration  of  the  definition  of  a  stable  system. 

The  curve  passing  through  P  is  tangent  to  L  and  represents 
the  orbit  of  a  galaxy. 
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tive  direction  on  the  line.  Further,  let  M  and  v  be  any  possible  values  of  the 
absolute  magnitude  and  of  the  peculiar  velocity  of  a  galaxy. 

Definition.  A  system  of  galaxies  will  be  called  stable  if,  whatever  be  P,  L,  M 
and  v,  the  'probability  that  a  galaxy  which  is  at  P  with  its  orbit  tangent  to  L  and  with 
its  absolute  magnitude  and  peculiar  velocity  equal,  respectively,  to  M  and  v,  will  be 
moving  in  the  positive  direction  of  L  equals  the  probability  that  it  will  be  moving  in 
the  negative  direction. 

Now,  using  formula  (5)  it  will  be  seen  that  in  a  stable  system  of  galaxies,  that 
is,  when  0  =  0  and  the  random  peculiar  velocity  satisfies  the  definition  just 
given,  for  every  given  P,  L,  M,  v  and  for  any  given  value  of  |cos  \p\,  the  frequency 
that  cos  ^  is  positive  will  be  equal  to  the  frequency  that  it  is  negative.  Then,  it 
follows  from  (5)  that  E(p  —  HD  cos  <p \p,  ip)  is  zero  so  that  for  a  given  the  radial 
velocity  p  is  not  correlated  with  p.  On  the  other  hand,  the  two  test  criteria  T 
and  Z  represent,  essentially,  the  least  squares  criteria  for  testing  the  hypothesis 
that  the  regression  coefficient  of  p  —  HD  cos  <p  on  m  (or,  equivalently,  on  p)  is 
equal  to  zero.  Thus,  if  the  systems  of  galaxies  considered  really  are  stable  in  the 
sense  of  the  above  definition  and  if  the  number  of  observations  is  sufficient  for 
the  least  squares  estimates  0  to  be  approximately  normally  distributed,  then  the 
criteria  T  and  Z  will  reject  the  hypothesis  tested  with  relative  frequencies 
approximately  equal  to  the  chosen  level  of  significance  a. 

From  the  point  of  view  of  the  frequency  with  which  the  two  tests  indicate 
instability  when  the  systems  of  galaxies  are  in  fact  stable,  we  see  that  the  two 
tests  are  approximately  valid  under  conditions  much  broader  than  those  en¬ 
visaged  in  their  deduction.  Unfortunately,  this  is  not  true  when,  in  fact,  the 
systems  are  unstable.  Should  it  happen  that,  for  example,  the  inner  part  of  a 
system  contracts  while  the  outer  part  expands  then,  as  is  easy  to  visualize,  the 
statistic  0  will  be  an  estimate  of  the  average  value  of  the  expansion  factors, 
appropriate  to  the  various  distances  from  the  center,  and  this  average  may 
happen  to  be  close  to  zero.  In  this  way,  the  two  tests  may  fail  to  detect  instability. 

6.  Asymptotic  power  of  the  two  tests  for  stability 

The  power  of  a  test  of  a  statistical  hypothesis  means  the  probability  that  this 
test  will  reject  the  hypothesis  when,  in  fact,  the  hypothesis  is  false.  Usually, 
there  are  many  different  ways  in  which  the  hypothesis  tested  may  be  false  and, 
naturally,  the  value  of  the  power  will  depend  on  the  kind  and  on  the  degree  of 
the  falsehood  of  the  hypothesis  under  test.  It  is  obvious  also  that,  in  order  to 
calculate  the  probability  representing  the  power  of  a  test,  it  is  necessary  to 
begin  with  rather  specific  conditions  for  the  problem.  In  this  section  we  return 
to  hypotheses  (i)  to  (vi)  adopted  in  section  2. 

The  hypothesis  of  stability  of  s  systems  of  galaxies  is  represented  by  s  equa¬ 
tions  0i  =  02  =  •  •  •  =  0s  =  0  indicating  that  the  postulated  expansion  coeffi¬ 
cient  in  each  of  the  systems  is  equal  to  zero.  This  hypothesis  may  be  false  in 
many  ways  including  the  following: 
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(a)  all  the  systems  expand  or  contract  at  exactly  the  same  rate  0  7^  0,  so  that 
0x  =  02  =  *  *  •  =  0,  =  0.  In  this  case  the  power  of  the  tests  considered  must  be  a 
function  of  the  common  rate  of  instability  0. 

(b)  The  rates  of  expansion  of  particular  systems  need  not  be  equal  but  at 
least  some  of  them  differ  from  zero.  In  this  case  the  power  of  the  tests  will  de¬ 
pend  upon  the  s  arguments  0i,  02,  •  •  •  ,0*. 

Case  (b)  is  more  general  than  (a)  and  therefore  in  our  study  we  shall  consider 
(b)  and  then  revert  to  (a)  as  a  particular  case. 

Reviewing  the  basic  hypotheses  (i)  to  (vi)  it  will  be  seen  that  they  are  not 
completely  specific  and  fail  to  indicate  the  orders  of  a  number  of  parameters 
involved.  The  parameters  Dt  and  c\  represent  one  category  of  “nuisance  param¬ 
eters”  and  the  two  tests  are  adjusted  so  that,  when  the  hypothesis  tested  is  true, 
it  will  be  rejected  with  a  preassigned  probability  a,  irrespective  of  the  true  values 
of  the  Di  and  ci.  However,  the  power  of  the  test  is  likely  to  depend  upon  the 
actual  values  of  these  parameters.  Further,  the  hypotheses  (i)  to  (vi)  visualize 
four  more  parameters,  namely  c,  cm,  t  and  H,  which  are  treated  as  known  num¬ 
bers,  determined  by  some  other  studies.  It  was  found  above  that  the  values  of 
the  two  criteria  are  independent  of  H  and  are  not  influenced  very  much  by 
moderate  errors  that  may  be  involved  in  the  assumed  value  of  c,  cm  and  r. 
However,  this  need  not  be  the  case  of  the  power.  In  fact,  the  power  is  likely  to 
depend  on  the  true  values  of  c,  cm,  t  and  H  and  also  on  those  assumed  values  of 
the  parameters  c,  cm  and  t  that  are  used  in  computing  the  test  criteria.  For  this 
reason  it  is  necessary  for  us  to  make  a  clear  distinction  between  the  two  sets  of 
values.  Thus,  the  original  symbols  a,  cm ,  and  r  will  be  used  to  denote  the  true 
values  of  the  dispersion  of  a  single  coordinate  of  the  position  of  a  galaxy  within 
a  system,  of  the  dispersion  of  the  absolute  magnitude  and  of  the  absorption 
coefficient.  These  true  values  are  unknown  but  this,  of  course,  should  not  pre¬ 
vent  us  from  discussing  them.  Along  with  the  true  parameter  values  c,  cm  and  r, 
we  shall  consider  their  “presumed”  values  used  in  the  application  of  the  tests. 
These  presumed  values  will  be  denoted  by  a,  cm  and  r,  respectively. 

It  is  known  (see,  for  example,  [17])  that  the  power  function  of  a  test  based  on 
a  criterion  of  the  type  T  may  be  approximately  calculated  from  the  formula 

1  f +»{a)-Q 

(47)  0(Q)  =  l-—==  e-*V2  dx, 

V2ir  J  —  v(<x)  —  q 

or 

(48)  A(Q)  -  4-  /“  e-^dx, 

V  27 T  Jv(2a)-Q 

depending  upon  whether  the  test  is  used  against  “two-sided”  or  “one-sided” 
alternatives,  where  v(a)  is  adjusted  to  the  adopted  level  of  significance  a  so  that 
0(0)  =  a,  and  Q  depends  on  the  extent  and  manner  in  which  the  hypothesis 
tested  is  false.  The  limitations  on  formulas  (47)  and  (48)  include  certain  con¬ 
ditions  of  regularity  and  the  assumptions  that  the  hypothesis  tested  is  only 
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“slightly”  false  while  the  number  of  observations  is  considerable.  It  is  to  this 
particular  assumption  that  the  description  “asymptotic”  refers. 

It  is  seen  that,  for  the  two-sided  alternatives,  assuming  0  <  0  and  0  >  0  as 
possibilities,  /3(Q)  is  an  increasing  function  of  |Qj.  For  the  one-sided  case  ft. (Q)  is 
increasing  with  Q.  It  follows  from  (47)  and  (48)  that  in  order  to  compute  the 
approximate  value  of  the  power  of  the  T-test  it  is  sufficient  to  compute  Q, 
which  must  depend  on  the  values  of  0i,  02,  •  •  •  ,  0„,  on  the  true  values  of  the 
nuisance  parameters  Z),  and  ai  and  on  the  true  values  of  the  remaining  four 
parameters,  a,  <tm,  t,  H  and  on  the  “presumed”  values  a,  vm,  and  f.  The  process 
of  computing  Q  consists  in  making  certain  substitutions  in  the  formula  repre¬ 
senting  the  test  criterion  T,  as  follows. 

The  expression  for  T  in  its  exact  form  (18)  depends  on  estimates  <ti  and  Z)„ 
the  latter  entering  the  formula  through  the  symbols  Ai3.  Our  first  step  in  com¬ 
puting  Q  is  to  replace  these  estimates  by  the  true  values  of  the  parameters. 

Next  we  notice  that  the  quant  ities  A  a  depend  also  upon  “presumed”  values 
of  <r,  ctm  and  r,  that  is,  on  a,  ctm  and  f.  These  presumed  values  are  left  in  the 
expressions  for  the  A  and  these  symbols  so  calculated  will  be  denoted  by  Z,y. 
Correspondingly,  if  we  deal  with  the  approximation  (37),  we  shall  write 

(49)  Atj  -  I =  .a~XF2— 2  ) 

Om  T  A((T 

where  the  bars  indicate  the  use  of  the  “presumed”  values  of  the  parameters 
involved.  [For  the  sake  of  brevity,  the  symbol  m  in  (49)  denotes  the  apparent 
magnitude  corrected  for  redshift.] 

Our  last  step  in  obtaining  Q  from  T  consists  in  replacing  each  pt;-  by  its  ex¬ 
pectation  evaluated  with  the  use  of  the  true  parameter  values  that  one  wants 
to  consider.  Then 

(50)  E (p o | th ij ,  (pij)  =  (D i  d-  diA  ij )  cos  <pij 
or,  with  the  approximation  (37)  and  putting  cos  <p,y  =  1, 

(51)  E{pij\mi},  <pij)  =  DiA-  0%  [ Ai .  +  (w,7  -  m<.)  J, 

where  the  absence  of  bars  indicates  “true”  parameter  values.  As  a  result  of  these 
substitutions  and  after  some  easy  algebra,  we  obtain 


(52) 

where 

(53) 


where  Wi  represents  the  value  of  Wi  corresponding  to  the  presumed  values  of  a, 
(tm  and  r,  and  where  0  is  a  weighted  mean  of  the  0;  with  weights  equal  to  the  prod¬ 
ucts  w%Wi.  Note  that  wq  is  simply  the  approximation  to  (1  A)[2Z(^v  —  A».)2]1/2. 
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Formula  (52)  indicates  that,  if  some  of  the  systems  of  galaxies  expand  and 
some  contract,  the  value  of  Q  may  well  be  zero  and  the  power  of  the  T-test  may 
be  nil.  Otherwise,  for  a  fixed  0,  the  quantity  Q  is  an  increasing  function  of  a, 
roughly  proportional  to  the  linear  dimensions  of  the  system  of  galaxies  and  a 
decreasing  function  of  the  dispersion  in  redshift.  Finally,  it  is  important  to  notice 
that  the  greater  the  variability  of  the  apparent  magnitudes  of  objects  within 
any  given  system,  the  greater  the  power  of  the  test.  If  for  a  given  system  all 
the  Ui  objects  used  in  the  test  had  exactly  the  same  apparent  magnitude,  then 
the  contribution  of  this  system  to  the  power  of  the  test  would  be  exactly  zero. 

The  last  factor  in  (52)  is  of  interest.  By  the  Schwarz  inequality,  it  cannot 
exceed  the  limit 

X)  WiWi 

(54)  ~  w2<)m’ 

i 

and  this  limit  is  reached  when  Wi  is  proportional  to  ie,-.  In  order  that  this  pro¬ 
portionality  hold  irrespective  of  the  distances  Di}  it  is  necessary  and  sufficient 
that  the  presumed  value  f  equal  the  actual  absorption  coefficient  t  and  that  the 
dispersions  a  and  <tm  be  underestimated  or  overestimated  by  a  and  &m  in  the 
same  proportion 


In  all  other  cases,  in  particular  when  the  absorption  coefficient  t  is  not  guessed 
correctly,  the  power  of  the  T-test  will  be  smaller  than  it  could  be. 

Considering  the  same  formula  (52)  from  a  slightly  different  point  of  view  we 
notice  that,  because  of  the  dependence  X*  =  r  -f-  a/Diy  whatever  be  the  presumed 
values  of  the  parameters,  the  power  of  the  T-test  increases  with  an  increase  in 
the  true  absorption  coefficient  r,  which  confirms  the  remarks  of  Dr.  Minkowski, 
made  at  the  discussion  of  this  paper. 

The  same  expression  for  X,  indicates  that,  other  things  being  equal,  the  close- 
by  systems  of  galaxies  contribute  to  the  detecting  power  of  the  test  much  more 
than  the  distant  systems,  with  large  Z>t.  Finally,  the  contribution  of  a  system  to 
the  power  depends  very  much  on  the  sums  ~  mi-)2>  that  is,  on  the  number 

of  objects  for  which  there  are  data  and  also  on  the  degree  to  which  the  apparent 
magnitudes  (corrected  for  redshift)  of  these  objects  diverge  from  each  other. 
Here  again,  because  of  the  presence  of  selection  of  objects  bright  enough  for 
convenient  observation,  the  more  distant  systems  are  at  a  disadvantage. 

Whatever  was  just  said  about  the  power  of  the  T-test,  applied  to  s  >  1  dif¬ 
ferent  systems  of  galaxies,  could  be  repeated  with  reference  to  the  application 
of  the  same  test  to  any  single  system  and  then  transferred  to  the  power  of  the 
Z-test  applied  to  s  systems.  Let  Qi  denote  the  value  of  Q  in  (52)  referring  to 
just  one  system  of  galaxies,  namely  the  ith,  so  that 

2 

Qi  =  6i-  Wi . 

<r  l 


(56) 
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The  asymptotic  power  of  the  Z-test  is  represented  by  the  power  of  the  F- test 
with  degrees  of  freedom  fx  =  s  and  /2  =  E  ni  —  2s,  and  with  noncentrality 
parameter 

(57)  E  Q2i  =  il  el  wt . 

x  Ol  i 

It  will  be  noticed  that,  if  all  the  s  systems  of  galaxies  expand  or  contract  at 
the  same  rate  6,  then  formula  (57)  reduces  to  the  square  of  Q  in  (52)  at  its 
maximum  value  attainable  only  when  the  presumed  value  r  =  r  and  equation 
(55)  is  satisfied.  Coupled  with  the  fact  that  the  noncentrality  parameter  (57) 
does  not  depend  upon  the  signs  of  the  particular  0f,  this  circumstance  may 
suggest  that  the  Z-test  has  a  definite  advantage  over  the  T-test  which,  for  its 
top  efficiency,  requires  the  sameness  of  the  signs  of  the  0*  and  a  lucky  guess  as 
to  the  parameters  cr,  <jm  and  t.  Unfortunately,  the  situation  is  not  so  simple 
because  for  a  fixed  value  of  the  noncentrality  parameter  (57)  the  power  of  the 
Z-test  declines  markedly  with  an  increase  in  s.  This  phenomenon  is  illustrated 
in  table  I. 

TABLE  I 

Decrease  in  Power  of  Z-test  when  the  Noncentrality  Parameter 
Is  Kept  Constant  and  the  Number  s  of  Systems  Is  Increased 


E  m  =  120,  =  0.05,  E  Qi2  =  8.00 


s 

1 

2 

3 

4 

5 

6 

13 

•80 

•71 

•65 

•59 

•55 

•51 

Table  I  indicates  that,  if  on  a  priori  grounds  the  possible  instability  of  several 
systems  of  galaxies  is  expected  to  be  characterized  by  coefficients  6X  all  of  the 
same  sign,  and  if  one  is  confident  about  the  presumed  values  of  the  three  pa¬ 
rameters  <r,  <tm  and  r,  then  the  T-test  may  have  a  definite  advantage  over  the 
Z-test. 

Before  proceeding  to  the  empirical  part  of  this  study  we  shall  give  another 
formula  referring  to  the  case  where  the  “presumed”  values  a,  cm  and  f  are  not 
equal  to  the  true  parameter  values  a,  <tm  and  r.  This  formula  represents  the 
expected  value  of  0„  the  estimate  of  the  expansion  coefficient  of  the  ith  system 
of  galaxies.  Using  (25),  the  approximations  (39)  and  denoting  by  bars  the  inter¬ 
vening  quantities  based  on  presumed  values  of  the  parameters,  we  have 

(58)  E0i)  =  _2X.  m ,.)2  £  (mu  -  mx.)E{Pi]\xpih  m,y) 

i 

_  gjf  +  XjU2  \j 
V2  (Tm  +  M  O’2  Xi 

Previous  calculations  indicated  that  the  first  two  factors  are  practically  inde¬ 
pendent  of  the  distance  of  the  ith  cluster.  On  the  other  hand,  depending  upon 
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the  presumed  value  r  of  r,  the  dependence  of  the  third  factor  upon  the  distance 
Di  may  well  be  very  strong.  In  particular,  this  will  be  the  case  if,  as  in  most 
of  our  calculations,  the  presumed  value  f  =  0.  Then 

(59)  Tr1  +  iD< 

and,  if  t  >  0,  the  absolute  value  of  the  extimate  |0,-|  will  tend  to  overestimate 
|<h|  by  a  factor  which  is  increasing  with  the  distance  D».  On  the  assumption 
that  di  =  02  =  •  •  •  =  0,  =  6, 

(60)  E0i)  *kO(l  +  ZD<), 
where  k  stands  for  a2d2M/a2a-M. 

7.  Coma  Cluster 

As  mentioned  at  the  outset,  since  the  first  application  of  the  kinematic  test 
in  1953,  the  number  of  presumed  members  of  Coma  Cluster  for  which  the 
apparent  magnitude  and  the  radial  velocities  are  available  increased  from  12 
to  50.  Table  II  gives  the  totality  of  these  data  including  the  new  observations 
of  radial  velocity  by  Mayall  and  the  new  observations  of  magnitude  by  G.  E. 
Kron  and  M.  F.  Walker.  We  are  very  grateful  to  Drs.  Kron  and  Walker  for 
allowing  us  to  list  and  to  use  their  unpublished  observations. 

At  the  end  of  table  II  are  fisted  five  galaxies  for  which  the  radial  velocity 
has  been  observed  but  no  observations  of  magnitude  are  available.  These  gal¬ 
axies  cannot  be  used  m  the  present  study. 

In  table  II  the  first  column  gives  the  NGC  and  IC  catalogue  numbers,  when 
it  was  possible  to  decide  with  some  assurance  that  the  identification  was  correct. 
Those  galaxies  designated  “Anon”  were  in  most  cases  given  arbitrary  numbers 
in  an  observing  fist  of  galaxies  for  which  it  was  judged  that  spectrograms,  suit¬ 
able  for  redshift  measurements,  could  be  obtained  with  exposures  of  \\  to  3 
hours  with  the  Crossley  nebular  spectrograph.  In  all  cases  Kodak  IlaO  emulsion 
was  used,  and  the  slit  width  was  0.1mm,  or  about  4  seconds  of  arc.  Except 
when  the  slit  was  placed  across  a  close  pair,  it  was  oriented  east-west  in  position 
angle,  and  was  kept  fixed  on  the  objects.  The  positions  in  the  second  and  third 
columns  resulted  from  measurement  of  a  48-inch  Palomar  schmidt  plate  of  the 
Coma  Cluster,  kindly  provided  by  Dr.  F.  Zwicky.  In  the  reductions  the  cata¬ 
logued  positions  of  the  brighter  objects  were  used  to  obtain  the  scale  and  ori¬ 
entation  factors;  for  some  catalogued  galaxies  the  fisted  quantities  represent 
corrected,  or  slightly  adjusted  values.  The  fourth  column  gives  the  radial  distance 
for  a  center  referred  to  the  midpoint  of  the  fine  joining  the  two  brightest  mem¬ 
bers,  NGC  4874  and  4889.  The  measurements  were  made  by  using  dividers  on 
a  negative  paper  print,  obtained  from  a  film  positive  contact-exposed  through 
the  glass  negative.  Although  these  distances  are  given  to  O'l,  the  larger  values 
may  be  systematically  in  error  by  greater  amounts,  since  only  a  linear  scale 


Redshifts  are  corrected  for  galactic  rotation  by  300  cos(i  —  55°)cos  b 
Magnitudes  are  not  corrected  for  latitude  effect  or  if-correction 
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nr.  IC  3960  56.7  28  08  10.9  +6880  +23828  17.4  [20] 

NGC  4865  56.8  28  21  9.7  +4655  -  3315  14.7  [13] 

NGC  4869  56.8  28  11  5.8  +6715  -  4657  (14.6)  [18]  (0.90) 
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NGC4676  12h44“5  +31°00'  244.1  +6564  (from  [22]) 

NGC  4858  56.7  28  22  13.8  +9398 

NGC  4922  12  59.1  29  35  83.4  +7372 

NGC  4944  13  01.5  28  28  54.4  +7013 

NGC  4961  13  03.5  28  00  80.4  +2579 
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factor  was  used.  The  new  redshifts  in  the  fifth  column  were  obtained  in  the 
same  way  as  those  reported  in  [13],  and  they  may  be  regarded  as  of  comparable 
accuracy;  a  colon  (:)  indicates  either  a  weak  plate  or  lines  of  poor  visibility  for 
measurement.  The  new  magnitudes  in  the  eighth  column,  identified  by  [18] 
in  the  ninth,  and  the  new  colors  in  the  last  column  resulted  from  two-color 
observations,  with  a  1P21  photomultiplier,  filters  that  define  the  P,  V  system, 
and  an  aperture  of  44  seconds  of  arc.  This  size  was  large  enough  to  include 
nearly  all  the  fight  for  the  majority  of  the  newly  observed  galaxies;  magnitudes 
in  parentheses  represent  results  from  single-night  observations.  The  value  of 
A  given  in  the  seventh  column  were  computed  from  formula  (30)  or  (35). 

Among  the  new  redshifts  there  are  four  that  seem  to  be  “abnormally”  low: 
NGC  4670,  Anon  1,  NGC  4725  and  4793.  Of  these,  NGC  4725  is  by  far  the 
largest  and  brightest,  and  in  fact  is  2  \  to  3  magnitudes  brighter  than  the  two 
brightest  galaxies  used  to  define  the  cluster  center.  It  is  also  more  than  3° 
distant  from  the  center  of  the  cluster.  Although  the  other  three  galaxies  have 
apparent  magnitudes  within  the  range  of  all  the  others,  their  colors  of  +0.32, 
+0.57 :  and  +0.48  mag  are  among  the  bluest  of  those  measured  in  this  region. 
Moreover,  the  bluest,  NGC  4670,  has  a  strong  emission  spectrum  and  is  of 
irregular  and  mottled  appearance,  while  both  Anon  1  and  NGC  4793  also  have 
emission  3727  [Oil],  together  with  an  early-type  continuum  in  the  UV.  These 
characteristics  all  suggest  dwarf  subluminous  systems.  Thus  these  four  galaxies 
of  notably  smaller  redshift  could  reasonably  be  considered  as  foreground  objects. 
On  the  other  hand,  these  newer  data  do  not  suggest  the  presence  of  any  back¬ 
ground  objects:  none  of  those  observed  have  redshifts  of  9500  km/sec  or  larger. 
This  result  is  qualitatively  to  be  expected,  since  in  a  given  solid  angle  the 
probability  of  encountering  the  more  numerous  nearby  dwarf  galaxies  is  greater 
than  that  of  finding  the  more  infrequent  giant  systems. 

The  first  panel  in  figure  3  gives  the  scatter  diagram  of  50  points  exhibiting 
the  relationships  between  the  radial  velocity  and  the  apparent  magnitude.  It 
will  be  seen  that  46  objects  provide  a  more  or  less  cohesive  group  of  points  in 
the  scatter  diagram,  while  the  remaining  four  represent  “outliers,”  already 
mentioned  at  the  beginning  of  the  paper.  Whether  these  outliers  should  be 
treated  as  members  of  the  Coma  Cluster  or  not  is  a  subjective  matter  and,  in 
order  to  allow  the  reader  to  follow  his  own  intuition,  the  kinematic  test  for 
stability  of  the  Coma  Cluster  has  been  applied  twice,  once  including  all  50 
objects  and  next  with  the  omission  of  the  outliers. 

Using  data  of  table  II  and  formulas  (41),  we  compute  the  sample  variance 
of  the  apparent  magnitude,  not  corrected  for  redshift,  its  covariance  with  the 
radial  velocity  and  the  sample  variance  of  the  latter.  Then  formulas  (43)  give 
the  sample  covariance  and  variance  of  the  apparent  magnitude  corrected  for 
redshift.  The  results  are  given  in  table  III. 

The  sample  variances  of  the  radial  velocities  represent  squares  of  1745  and 
932  km/sec,  respectively.  It  should  be  noted  that  they  bracket  the  value  of 
1000  km/sec  estimated  by  Oort  [10]  using  23  galaxies  in  the  central  region  of 
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TABLE  III 
Coma  Cluster 

Variances  and  Covariances  of  Radial  Velocity  and  Apparent  Magnitude 


n 

p. 

S2(m) 

S(m,  p) 

S2(p) 

o r{m,  p) 

a2(m) 

With  outliers 
Without 

50 

6432 

1.000 

682.27 

3,044,423 

636.61 

0.980 

outliers 

46 

6866 

0.481 

-102.07 

868,448 

-115.10 

0.484 

the  cluster.  Also,  working  backwards  from  the  assumption  of  stability,  Oort 
found  the  value  583  km/sec. 

In  order  to  proceed  further,  we  must  now  adopt  some  “presumed”  values  for 
the  three  parameters  <r,  <tm,  and  r  for  which  no  independent  estimates  are  at¬ 
tempted  in  the  present  study.  The  presumed  value  of  cr  used  below  is  = 
6.866  X  10 5/h.  It  was  obtained  by  studying  Shane’s  counts  of  presumed  members 
of  the  Coma  Cluster  made  in  rings  about  the  apparent  cluster  center.  The  dis¬ 
tribution  of  these  counts  represents,  approximately,  the  distribution  of  sin  <p, 
with  =  a/D.  The  value  of  found  is  of  the  order  of  0.01.  Using  the  average 
radial  velocity  of  presumed  members  of  Coma  Cluster  (without  outliers)  given 
in  table  III  and  the  estimate  of  distance  D  =  10 *p-/h,  and  solving  for  <r  the 
above  value  for  a  is  obtained.  Subsequent  studies  of  Shane  [19]  indicate  that 
a  —  6.9  X  10 6/h  is  reasonably  consistent  with  his  counts  relating  to  other 
clusters. 

The  adopted  “presumed”  value  of  <tm  is  om  -  1.25  mag,  which  is  a  compro¬ 
mise  of  the  various  current  estimates.  The  “presumed”  value  of  r  used  is  f  =  0. 

Using  these  values  further  computations  yield  the  results  exhibited  in  table  IV. 


TABLE  IV 

Kinematic  Test  for  Stability  of  the  Coma  Cluster 


Conditional 
redshift 
dispersion 
<ti  (km/sec) 

Expansion 

coefficient 

8(km/sec/pc) 

S.E.  of  exp. 
coefficient 

Test 

criterion 

T  —  0/a$ 

P\T) 

With  outliers 
Without 

1639 

+0.0636  h 

0.0232  h 

+2.740 

0.009 

outliers 

927 

-0.0250  h 

0.0208  h 

-1.201 

0.236 

The  last  column  of  table  IV  gives  the  probability  that,  under  conditions  of 
stability  0  =  0  and  due  to  chance  alone,  the  absolute  value  of  the  criterion  T 
will  equal  or  exceed  the  value  actually  obtained  for  the  Coma  Cluster,  as  given 
in  the  preceding  column.  Also,  this  is  the  upper  bound  of  the  probability  of 
being  wrong  if  we  make  a  rule  of  rejecting  the  hypothesis  tested  in  conditions 
similar  to  those  referring  to  the  Coma  Cluster.  It  is  seen  that  if  the  four 
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outliers  are  bona  fide  members  of  the  Coma  Cluster,  the  risk  in  asserting  that 
the  cluster  is  unstable  is  very  small.  The  positive  sign  of  the  estimate  of  0  indi¬ 
cates  expansion.  However,  if  the  outliers  are  excluded,  which  the  present  authors 
think  appropriate,  then  the  picture  changes  entirely.  First,  the  upper  bound  of 
the  probability  of  error  in  asserting  instability  increases  to  almost  one  in  four. 
Second,  the  estimate  of  0  obtained  with  outliers  omitted  is  negative  and  indicates 
contraction. 

In  judging  these  results  it  must  be  borne  in  mind  that  the  limit  (the  so-called 
level  of  significance  a)  between  “small”  and  “large”  probabilities  of  error  is 
conventional;  it  is  a  subjective  matter  whether  to  risk  an  assertion  when  the 
probability  of  error  may  be  as  large  as  24  per  cent.  The  second  relevant  con¬ 
sideration  is  that,  if  the  level  of  significance  is  set  low,  the  chance  of  detecting 
an  effect,  when  it  exists  and  is  reasonably  large  (that  is,  the  power  of  the  test) 
is  small.  In  cases  where  observational  data  are  not  very  difficult  to  obtain,  the 
usual  values  of  the  level  of  significance  are  0.05  and  0.01.  With  the  high  variabil¬ 
ity  of  the  radial  velocity,  in  order  to  have  a  reasonable  chance  of  detecting  ex¬ 
pansion  working  with  a  =  0.05,  or  even  a  =  0.10,  the  necessary  amount  of 
data  would  be  prohibitive  (see  table  X).  For  these  reasons,  the  present  authors 
are  inclined  to  adopt  a  much  more  liberal  level  of  significance  and  to  consider 
that  values  P{T}  ^  0.24  are  a  slight  indication  of  instability  to  be  taken  into 
account  when  considering  other  evidence.  However,  one  should  not  be  very 
much  surprised  if  the  present  indication  of  contraction  of  the  Coma  Cluster  is 
contradicted  by  further  studies.  Our  final  conclusion  is  that,  if  the  four  outliers 
are  members  of  the  Coma  Cluster,  this  cluster  must  be  expanding.  If  not,  then 
there  is  slight  evidence  of  contraction. 

With  a  casual  approach  to  the  problem  of  stability  or  instability  of  systems 
of  galaxies,  one  is  inclined  to  treat  the  phenomenon  in  its  extreme  simplicity :  a 
cluster  or  a  group  may  be  stable  or,  if  not,  it  may  either  expand  as  a  whole  or 
contract  as  a  whole.  However,  the  phenomenon  of  instability  of  systems  of  galax¬ 
ies  need  not  be  simple.  In  fact,  as  suggested  by  Oort  [10]  and  subsequently  re¬ 
affirmed  by  van  Albada  [14],  in  certain  phases  of  evolution  of  a  cluster  its 
inner  parts  may  be  contracting  and  its  outer  parts  may  be  expanding.  It  is, 
then,  not  impossible  that  the  “outliers”  are  actual  members  of  the  Coma  Cluster 
but,  being  located  in  its  outer  regions,  are  in  the  process  of  “evaporating.”  In 
this  case  the  quantities  0  obtained  from  the  data  are  estimates  of  average  rates  of 
instability,  one  relating  to  the  cluster  including  the  outer  regions,  and  the  other 
restricted  to  the  inner  core.  On  this  particular  assumption,  the  velocities  of 
“evaporation”  of  cluster  members  in  the  outer  regions  must  be  substantially 
higher  than  might  be  computed  using  the  estimate  given  in  the  table. 

With  observations  of  only  50  galaxies  available,  any  subdivision  of  the  data 
according  to  the  value  of  will  give  shaky  results.  Nevertheless,  it  seemed 
interesting  to  carry  out  the  test  when  the  observations  were  divided  into  three 
rings:  (i)  <  0?2  with  12  galaxies,  (ii)  0!2  <  <p  <  0!5  with  14  galaxies,  and 

(iii)  (p  >  0?5  with  the  remaining  24  galaxies.  We  find  that  for  the  two  inner 
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rings  there  is  no  evidence  of  instability.  This  is  not  surprising  with  the  small 
value  of  n.  On  the  other  hand,  the  results  for  the  outermost  ring  almost  duplicate 
those  for  the  entire  cluster,  with  strong  evidence  for  expansion  or  mild  for 
contraction  according  to  whether  the  outliers  are  included  or  not. 

The  above  discussion  emphasizes  the  difficulties  which  are  connected  with 
efforts  to  study  individual  systems  of  galaxies  when  a  substantial  effort  is  made 
to  collect  data  relating  to  members  of  a  particular  system.  Here  a  subjective 
judgment  is  unavoidable  and  may  easily  lead  to  spurious  conclusions.  The  alter¬ 
native  approach  is  to  give  up  entirely  the  idea  of  identifying  members  of  par¬ 
ticular  systems  and  to  conduct  the  study  on  data  for  all  galaxies  in  a  given 
region  of  the  sky  selected  on  some  well  defined  basis,  for  example,  on  the  basis 
of  estimated  ease  in  measuring  radial  velocity.  Unfortunately,  studies  of  this 
kind  require  a  much  more  complicated  theory  with  more  parameters  to  estimate 
and,  therefore,  considerably  greater  volume  of  observational  data. 

The  actual  values  of  the  two  estimates  0  given  in  table  IV  indicate  that  the 
presumed  values  a,  &m  and  t  =  0  are  wrong  causing  the  coefficient  of  0  in  for¬ 
mula  (60)  to  be  much  greater  than  unity.  In  fact,  if  the  actual  value  of  |0|  were 
of  the  order  of  0.06,  then  the  “velocity  of  instability,”  that  is,  of  expansion  or 
of  contraction,  of  a  galaxy  at  a  typical  distance  from  the  center  of  the  cluster 
equal  to  <r  =  6.9  X  105  pc  would  have  been  of  the  order  of  45,000  km/sec, 
which  is  fantastic.  Although  the  subject  of  the  present  paper  is  limited  to  testing 
the  hypothesis  that  0  =  0  and  the  problem  of  estimating  the  value  of  0  is  not 
our  immediate  concern,  it  is  of  some  interest  to  consider  the  effect  on  the  estimate 
of  0  of  some  plausible  errors  in  the  presumed  a,  aM  and  f . 

The  adopted  a  is  very  poorly  determined  and  we  should  not  be  surprised  to 
find  that  the  actual  a  =  2a.  Similarly,  as  indicated  in  [16],  aM  may  well  be  close 
to  unity.  Finally,  the  studies  of  Zwicky  suggest  that  t  may  be  of  the  order  of 
10-7  mag/pc.  Taking  r  =  2  X  10~7  mag/pc  we  find 


(61) 


of  +  Vo2  / 1  I  D 
a2  +  XV  \  a 


13. 7\ 
2.17,/ 


41.5. 


The  combination  of  the  above  errors  would  then  reduce  the  estimates  of  0  by  a 
factor  of  about  40,  yielding  a  value  only  about  15  times  larger  than  the  Hubble 
constant.  At  the  distance  a  from  the  center  of  the  cluster  this  would  imply  the 
velocity  of  instability  of  the  order  of  1000  km/sec,  still  very  large  but  not  as 
exorbitant  as  on  the  original  presumed  values  of  the  parameters. 

The  necessity  of  revising  the  originally  adopted  presumed  values  of  a,  9m  and 
f  =  0  is  also  indicated  by  the  study  of  the  power  of  the  T-test.  Table  V  gives 
the  approximate  value  of  the  probability  that  the  T-test  applied  to  the  data  for 
the  Goma  Cluster,  without  the  outliers,  will  indicate  significance  at  24  per  cent 
(the  observed  value  of  the  probability  as  shown  in  table  IV),  computed  on 
several  hypotheses  regarding  the  intervening  parameters. 

It  is  seen  that,  if  the  presumed  values  9m  and  a  were  entirely  correct  and  the 
expansion  coefficient  0  coincided  with  the  Hubble  constant,  then,  even  with 
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TABLE  V 

Asymptotic  Power  of  7-test  Applied  to  Data  on  Coma  Cluster, 

Without  Outliers 
In  all  cases  it  is  assumed  <tm  =  9m 


Hypothetical 
|0|  and  a 

T  =  0 

Hypothetical  r 
r  =  10~7 

r  =  2  X  10~7 

l«l  =H 

0.24 

0.24 

0.24 

a  =  9 

|0|  =10  H 

0.25 

0.42 

0.68 

a  =  29 

t  —  2  X  10~7,  the  chance  of  obtaining  the  observed  T  would  be  as  small  as 
with  6  =  0.  The  assumption  that  |0|  equals  ten  times  the  Hubble  constant, 
combined  with  the  assumption  that  the  actual  value  of  a  is  twice  the  “pre¬ 
sumed"  and  with  the  assumption  t  =  0,  does  not  help  much.  It  is  only  when 
we  assume  that  r  =  10~7  combined  with  |0|  =  10  H  and  <r  =  2 a  that  we  reach 
nonnegligible  chance  that  the  test  will  detect  instability.  Since,  in  fact,  indica¬ 
tions  of  instability  have  been  detected,  our  conclusion  is  that,  unless  the  data 
available  are  dominated  by  effects  of  selection  distorting  the  distribution  of  the 
radial  velocity,  the  figures  in  table  IV  are  indicative  not  only  of  instability 
but  also  of  the  presence  of  a  substantial  amount  of  absorbing  matter  within 
the  Coma  Cluster. 

8.  Joint  study  of  seven  clusters  and  fifteen  groups  of  galaxies 

In  the  present  section  we  give  the  results  of  the  application  of  the  kinematic 
test  for  stability  to  all  clusters  and  all  groups  of  galaxies  for  which  we  could 
secure  unambiguous  data  on  apparent  magnitudes  and  radial  velocity  for  at 
least  three  member  objects.  Most  of  the  data  come  from  the  fundamental  paper 
by  Humason,  Mayall  and  Sandage  [13].  Here  we  believe  to  have  exhausted  the 
information  given  in  their  table  II  “Redshift  from  Bright  Nebulae  in  Clusters" 
for  which  we  could  find  apparent  magnitudes  in  table  Al.  We  had  somewhat  less 
luck  with  the  data  of  table  III  “Redshifts  from  Faint  Nebulae  in  Clusters," 
where  quite  a  few  clusters  had  to  be  omitted.  One  example  is  the  Hydra  Cluster 
for  which  three  radial  velocities  are  given.  However,  one  of  these  radial  velocities 
is  from  a  blend  of  spectra  of  two  different  galaxies  and  it  was  judged  appropriate 
to  discard  it.  Also  in  a  number  of  cases,  we  found  it  difficult  to  match  the  radial 
velocities  with  apparent  magnitudes.  As  a  result  of  all  such  difficulties,  out  of 
the  18  clusters  treated  in  [13]  we  could  use  only  six.  The  radial  velocities  for  the 
seventh  cluster  “around  NGC  6166"  were  communicated  to  us  by  Dr.  R. 
Minkowski  and  the  apparent  magnitudes  by  Dr.  G.  O.  Abell,  in  both  cases  in 
advance  of  publication,  and  it  is  a  pleasure  to  record  here  our  indebtedness. 
Also,  Dr.  Abell  was  kind  enough  to  supply  us  with  unpublished  apparent  mag- 
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nitudes  of  several  objects  in  the  Hercules  and  in  the  Corona  Borealis  Clusters. 
Finally,  we  wish  to  thank  Dr.  E.  M.  Burbidge  for  her  kind  help  in  matching  the 
apparent  magnitudes  and  the  radial  velocities  for  several  objects  in  the  Hercules 
Cluster. 

Most  of  the  data  for  groups  of  galaxies  stem  from  table  XI  “Data  for  Repre¬ 
sentative  Groups  of  Nebulae”  [13],  a  most  convenient  table  giving  side  by  side 
the  radial  velocities  and  the  apparent  magnitudes  for  each  galaxy.  These  data 
were  supplemented  by  those  for  Messier  81  and  Messier  101  as  published  by 
Holmberg  [23].  The  five  panels  of  figure  3  give  the  scatter  diagrams  of  radial 
velocities  and  apparent  magnitudes  for  Coma,  Virgo,  Ursa  Major  Cloud,  Her¬ 
cules,  and  G383.  It  will  be  noticed  that  the  Ursa  Major  Cloud  has  one  “outlier.” 
While  the  radial  velocity  of  the  particular  object  is  comparable  to  that  of  other 
presumed  members  of  this  group,  its  apparent  magnitude  is  much  brighter. 

The  basic  statistics  are  given  separately  for  clusters  and  for  groups  in  tables 
VI  and  VII,  respectively.  In  each  case,  the  systems  of  galaxies  are  arranged  in 


TABLE  VI 

Basic  Statistics  for  Seven  Clusters  of  Galaxies 
Magnitudes  corrected  for  latitude  effect  but  not  for  A’-correction 


Cluster 

n 

P. 

km/sec 

m. 

S(m) 

S(m,  p) 

S(p) 

Virgo  whole 

80 

1197 

11.4 

92.4 

641 

Virgo  without  S.E. 

70 

11.5 

102.9 

676 

Perseus 

4 

5516 

13.2 

Coma  whole 

50 

6432 

14.7 

682.3 

1745 

Coma  without  outliers 

46 

6866 

14.8 

-102.1 

932 

Around  NGC  6166 

19 

9028 

16.3 

93.7 

864 

Hercules 

15 

10776 

15.2 

0.52 

40.2 

652 

Pegasus  II 

3 

12821 

15.0 

-103.0 

662 

Cor.  Bor. 

8 

21651 

17.4 

-  84.2 

1294 

the  order  of  increasing  mean  corrected  radial  velocities.  The  first  column  gives 
the  number  n  of  objects  for  which  the  necessary  data  are  available  and  here  we 
must  express  our  regret  that  for  many  systems  this  number  is  pitifully  small. 
The  next  two  columns  give  the  mean  radial  velocity  and  the  mean  apparent 
magnitude  exhibiting  the  familiar  redshift-magnitude  relation.  The  next  column 
gives  the  sample  dispersion  of  the  apparent  magnitude  not  corrected  for  redshift. 
This  column  is  interesting  because  of  the  marked  negative  correlation  between 
redshift  and  S(m)  which  clearly  indicates  the  effect  of  selection  of  objects  bright 
enough  for  convenient  measurement  of  radial  velocities.  As  mentioned  earlier, 
this  selection  does  not  invalidate  the  kinematic  tests  for  stability.  However,  as 
was  pointed  out  before,  the  smaller  the  variability  of  the  apparent  magnitude  of 
objects  of  a  given  system,  the  smaller  is  the  contribution  of  this  system  to  the 
power  of  the  tests.  The  next  column,  of  sample  covariances  of  the  radial  velocity 
with  the  apparent  magnitude  (not  corrected  for  redshift),  shows  a  regrettable 
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TABLE  VII 

Basic  Statistics  for  Fifteen  Groups  of  Galaxies 
Magnitudes  corrected  for  latitude  effect  but  not  for  ^-correction 


Group 

n 

P. 

km /sec 

m. 

S(m) 

S(m,  P ) 

S(p) 

M81 

7 

172 

10.0 

1.39 

-  14.0 

132 

M101 

6 

478 

10.5 

1.59 

-  29.5 

102 

G1023 

5 

513 

10.5 

0.87 

-  51.1 

257 

Leo  Group 

18 

788 

11.1 

0.98 

36.9 

264 

Ursa  Major  Cloud 

27 

1181 

11.4 

0.74 

146.4 

418 

U.M.  without  outliers 

26 

1208 

11.5 

0.56 

81.1 

403 

G3190 

3 

1222 

12.2 

0.58 

-  40.8 

70 

G1068 

3 

1604 

10.8 

1.05 

472.1 

453 

Fornax 

3 

1631 

11.0 

0.10 

-  20.2 

299 

G5846 

8 

1808 

12.2 

0.87 

69.2 

384 

G7619 

6 

3836 

13.3 

0.98 

102.8 

285 

G6928 

3 

4643 

13.8 

1.02 

-153.7 

308 

G507 

3 

4665 

13.1 

0.74 

-265.7 

416 

G383 

8 

5264 

14.4 

0.90 

334.9 

570 

Stephan’s  Quintet 

4 

6735 

14.2 

0.68 

-  50.3 

470 

G68 

5 

6785 

14.9 

0.66 

117.5 

458 

lack  of  regularity  including  the  signs.  On  the  other  hand,  some  regularity  is 
noticeable  in  table  VII  in  the  last  column  of  dispersion  in  redshift  for  groups  of 
galaxies.  Because  of  the  very  unequal  and  occasionally  very  small  number  of 
observations  on  which  these  dispersions  are  based,  in  order  to  see  the  tendency 
more  clearly,  it  is  advisable  to  combine  several  groups,  say  at  intervals  of  500 
km/sec  of  their  radial  velocities,  and  to  compute  a  “pooled”  measure  of  dis¬ 
persion  for  each  combination.  If  one  does  this,  one  obtains  an  unmistakable 
picture  of  a  systematic  increase  of  dispersion  roughly  up  to  the  mean  radial 
velocity  of  1000  km/sec.  Since  the  validity  of  the  tests  introduced  here  is  most 
sensitive  to  distortions  of  the  conditional  distribution  of  the  radial  velocity,  the 
tendency  indicated  deserves  attention,  particularly  because  of  the  following 
comments  by  Humason  taken  from  [13]  page  101 :  “Among  the  nebulae  desig¬ 
nated  ‘anonymous’  are  several  originally  observed  as  possible  members  of 
clusters,  but  later  found  to  be  field  nebulae  with  redshifts  not  in  agreement  with 
those  obtained  from  cluster  members.”  In  connection  with  relatively  large 
values  of  sample  dispersion  of  magnitudes  and  relatively  small  dispersions  of 
radial  velocities  indicated  for  the  four  presumably  nearest  groups,  one  is  led  to 
suspect  that  the  assignment  of  objects  to  these  groups  may  have  been  affected 
by  considerations  of  similarity  of  redshift,  so  that  the  distribution  of  this  char¬ 
acteristic  in  these  groups  became  distorted.  If  this  actually  happened,  then  the 
application  of  the  test  would  reflect  not  so  much  the  actual  state  of  affairs 
within  the  groups  concerned  as  the  arbitrary  distortions  made  in  the  process  of 
collecting  the  data.  The  five  groups  for  which  there  is  the  suspicion  of  such  dis- 
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<  ort  ions,  because  of  the  excessively  small  dispersion  of  redshift,  arc  M81,  M101, 
G1023,  Leo  Group  and  G3190. 

Because  of  these  considerations,  the  summary  treatment  of  groups  of  galaxies 
has  been  performed  twice,  once  taking  into  account  all  the  fifteen  groups  for 
which  the  data  are  available,  and  next  for  the  ten  groups  for  which  no  indication 
of  selection  for  redshift  is  noticeable. 

A  glance  at  the  last  columns  of  tables  VI  and  VII  is  sufficient  to  indicate  the 
desirability  of  treating  the  clusters  and  the  groups  of  galaxies  separately;  there 
is  little  doubt  that  the  variability  of  the  radial  velocity  in  the  systems  of  these 
two  broad  classes  is  not  the  same. 

Out  of  all  the  clusters  and  groups,  the  number  of  observations  available 
justifies  individual  treatment  only  for  Coma,  given  above,  and  for  the  Virgo 
Cluster.  Table  VIII,  analogous  to  table  IV,  gives  the  results,  separately  for  the 


TABLE  VIII 

Kinematic  Test  eor  Stability  of  the  Virgo  Cluster 


Conditional 

redshif  t 
dispersion 
<n(km/se(0 

Expansion 

coefficient 

0(km/sec/'pe) 

S.E.  of  exp. 
coefficient 

*8 

Test 

criterion 

T  =  0/aS 

P\T\ 

Whole  cluster 

639 

0.00156  h 

0.00131  h 

1.193 

0.237 

Without  S.E. 

674 

0.001 74  h 

0.00148// 

1.172 

0.246 

whole  cluster  and  separately  for  its  main  part,  excluding  the  Southern  Extension. 

The  sign  of  the  expansion  coefficient  is  positive  for  both  the  cluster  as  a  whole 
and  for  its  part  after  the  omission  of  the  Southern  Extension.  Thus,  if  there  is 
instability,  we  would  expect  expansion.  If  we  admit  a  priori  that  the  cluster  may 
be  expanding  or  may  be  contracting,  the  risk  of  asserting  instability  of  the  Virgo 
Cluster  is  just  under  one  in  four.  On  the  other  hand,  if  a  priori  we  exclude  the 
possibility  of  contraction,  then  the  risk  is  halved.  The  indication  of  instability 
of  the  Virgo  Cluster  equals  in  intensity  that  of  the  Coma  Cluster  without  out¬ 
liers.  However,  in  one  case  the  indicated  instability  is  expansion  and  in  the 
other  contraction.  In  both  cases  the  conclusion  depends  on  the  hypothesis  that 
the  collection  of  data  did  not  distort  the  conditional  distribution  of  radial 
velocity. 

Remark:  out  of  the  80  objects  available  for  the  Virgo  Cluster  ten  objects  are 
within  the  “Southern  Extension”  of  the  cluster.  It  seemed  advisable  to  perform 
the  analysis  twice,  first  treating  all  the  80  objects  as  one  cluster  and  then  re¬ 
stricting  the  calculations  to  those  70  objects  not  in  the  Southern  Extension. 
However,  the  inclusion  or  the  exclusion  of  the  Southern  Extension  does  not 
make  much  difference.  This  is  seen  in  tables  VI  and  VIII  and  the  agreement 
persists  in  further  calculations.  For  this  reason,  in  describing  further  results, 
only  the  figures  relating  to  the  Virgo  Cluster  as  a  whole  are  given. 
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The  joint  treatment  of  all  clusters  together  and  then  of  all  groups  together  is 
based  on  the  presumption  that  there  is  a  certain  amount  of  similarity  in  the 
systems  combined  into  a  single  class.  The  particular  aspect  of  similarity  is  the 
equality  of  the  residual  variance  of  the  radial  velocity.  Taking  the  available 
data,  so  to  speak  at  their  face  value,  this  presumption  is  far  from  confirmed.  In 
fact,  for  example,  the  data  for  the  Coma  and  for  the  Virgo  Clusters  taken  as  a 
whole,  or  with  the  indicated  omissions,  indicate  that  members  of  the  Coma 
Cluster  have  a  greater  dispersion  of  radial  velocity  than  those  of  the  Virgo 
Cluster.  The  motivation  for  the  joint  treatment  of  all  systems  belonging  to  each 
of  the  two  classes  is  the  suspicion  that  each  of  the  systems  must  have  been  sub¬ 
jected  to  some  selection  and  that  probably  the  intensity  of  this  selection  varied 
from  case  to  case.  In  these  circumstances  it  appears  interesting  to  see  the  results, 
so  to  speak,  averaged  over  the  whole  class  of  systems  of  galaxies.  In  addition, 
the  smallness  of  the  number  of  observations  for  particular  groups  makes  any 
analysis  other  than  by  a  summary  test  a  rather  hopeless  undertaking. 

Table  IX  gives  pooled  estimates  of  the  conditional  standard  deviation  ai  of 


TABLE  IX 

Pooled  Estimates  of  Conditional  Dispersion  of  Radial  Velocities 


Category  of 
Systems  of 
Galaxies 

Rate  of  Instability  Assumed 

Common  to  all  systems 

Possibly  Different  in  each 

Outliers 

included 

excluded 

included 

excluded 

ai  d.f. 

<7,  d.f. 

<7!  d.f. 

<7,  d.f. 

7  clusters 

15  groups 

10  groups 

1114.5  171 
366.2  93 

788.6  167 

361.2  92 

403.5  58 

1078.7  165 
346.4  79 

791.4  161 

348.6  78 

401.7  49 

the  radial  velocity  calculated  separately  for  clusters  and  groups,  with  and  with¬ 
out  outliers  and  under  the  two  different  assumptions,  one  that  the  rate  of 
expansion  in  all  systems  of  a  given  category  is  the  same  and  the  other  that  each 
system  of  galaxies  may  expand  or  contract  at  its  own  rate. 

It  is  seen  that,  while  the  omission  of  the  four  outlier  galaxies  in  the  Coma 
Cluster  considerably  influences  the  pooled  estimates  of  the  conditional  disper¬ 
sion  in  all  clusters,  the  removal  of  the  single  outlier  galaxy  in  the  Ursa  Major 
Cloud  produces  hardly  any  difference  in  the  pooled  estimate  of  <ti.  On  the  other 
hand,  of  course,  the  elimination  of  five  groups  suspected  of  being  affected  by 
selection  for  similarity  in  radial  velocities,  marked  with  unusually  small  sample 
dispersions  in  redshift,  substantially  increases  the  pooled  estimates  of  condi¬ 
tional  dispersion  in  p. 

The  computations  also  were  performed  on  all  groups  except  the  three  with 
very  small  values  of  S(p),  the  dispersion  in  radial  velocity.  The  results  are 
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always  intermediate  between  those  for  all  15  groups  and  for  10  groups.  They 
are  not  shown  in  table  IX  nor  in  any  later  table  to  save  space. 

Tables  X  and  XI  give  the  details  for  the  application  of  the  summary  T-  and 
Z-tests,  separately  for  clusters  and  for  groups.  The  conclusions  suggested  by 
table  X  are  as  follows : 

(i)  If  we  take  it  for  granted  that  the  phenomenon  of  instability  is  the  same 
for  all  clusters  and  is  measured  by  a  common  value  of  the  expansion  coefficient  0, 
then  there  is  a  slight  indication  of  expansion  if  outliers  are  included,  but  none 
when  they  are  excluded. 

(ii)  If  we  admit  the  possibility  that  particular  clusters  may  be  unstable, 
each  with  a  different  value  of  0,  then,  with  outliers  included,  there  is  strong 
evidence  that  at  least  some  of  them  (and  this  would  be  Coma)  are  unstable. 
In  this  case,  the  chance  of  error  in  asserting  instability  of  at  least  some  of  the 
clusters  is  less  than  one  in  a  hundred.  If  the  outliers  are  excluded,  no  such 
evidence  exists;  the  chance  of  error  in  asserting  instability  may  be  as  large  as 
three  in  four. 

(iii)  Naturally,  lack  of  clear-cut  evidence  of  instability,  does  not  amount  to 
evidence  of  stability.  The  power  of  the  test,  as  shown  in  the  last  columns,  is  low 
when  0  =  H  even  with  a  =  0.10.  In  this  connection  our  attention  is  attracted 
by  a  substantial  correlation  between  the  mean  radial  velocity  p.  and  the  absolute 
value  of  the  estimate  0  of  the  expansion  factor.  With  reference  to  formula  (60) 
this  correlation  suggests  that  some  clusters  may  be  expanding,  0»  >  0,  and  some 
others  contracting,  0*  <  0,  but  all  contain  absorbing  matter  explaining  the 
correlation  observed. 

Asterisks  in  table  XI  mark  the  five  groups  of  galaxies  for  which  table  VII 
indicated  the  suspicion  that  objects  assigned  to  them  were  those  with  similar 
radial  velocities,  so  that  the  variability  of  p  was  artificially  decreased.  It  is  for 
this  reason  that  the  computations  were  performed  more  than  once,  first  for  all 
groups,  next  for  the  ten  groups  not  marked  with  the  asterisks,  and  then  for  the 
twelve  groups  for  which  the  variance  in  p  was  not  too  extreme  (results  not 
shown) . 

Among  the  group  there  is  just  one,  the  Ursa  Major  Cloud,  with  any  “outlier” 
member.  The  figures  in  table  XI  refer  to  calculations  made  with  the  omission 
of  the  outlier. 

The  conclusions  suggested  by  table  XI  may  be  summarized  as  follows. 

(i)  As  far  as  the  Z-test  is  concerned,  its  results  are  very  much  independent 
of  whether  we  include  all  15  groups  or  only  10.  In  both  cases,  the  probability  of 
error  in  asserting  instability  may  be  as  high  as  one  in  six.  However,  it  will  be 
remembered  that  with  a  substantial  number  of  groups  included  and  with  a  very 
small  number  of  observations  per  group,  the  power  of  the  Z-test  is  low. 

(ii)  If  we  adopt  the  a  priori  attitude  that  for  all  the  groups  treated  the  insta¬ 
bility  parameter  0  must  have  the  same  value,  then  the  conclusions  suggested  by 
the  T-test  differ  sharply  on  whether  it  is  applied  to  all  the  15  groups  or  only  to 
those  without  the  asterisk.  In  the  first  case,  there  is  no  evidence  of  instability. 
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Asterisks  indicate  groups  suspected  of  being  affected  by  excessive  selection 
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In  the  second,  the  probability  of  error  in  asserting  instability  is  less  than  one  in 
fifty.  On  the  other  hand,  there  is  an  obvious  argument  against  giving  credence 
to  this  result.  By  omitting  from  the  data  those  items  about  which  the  investi¬ 
gator  has  an  intuitive  suspicion,  one  is  very  much  at  the  mercy  of  subjective 
judgment. 

(iii)  The  absolute  values  of  the  estimates  0,-  exhibited  in  table  XI  show  the 
same  kind  of  correlation  with  the  mean  radial  velocity  p.  as  in  table  X.  This  is 
an  interesting  phenomenon,  suggesting  both  instability  and  the  presence  of 
absorbing  matter  within  the  systems  of  galaxies. 


9.  Summary  and  conclusions 

(i)  As  understood  in  this  paper,  the  stability  of  a  system  of  galaxies,  a 
cluster  or  a  group,  means,  roughly  speaking,  that  for  any  given  distance  from 
the  center  C  of  the  system  and  for  any  given  peculiar  velocity  v,  the  galaxies 
approach  C  with  this  velocity  as  frequently  as  they  move  away  from  it.  This 
definition  implies  that  for  members  of  a  stable  system  there  shall  be  no  correla¬ 
tion  between  the  redshift  and  the  apparent  magnitude,  corrected  for  redshift. 

(ii)  Two  locally  and  asymptotically  most  powerful  tests  for  stability  of 
several  systems  of  galaxies  have  been  introduced,  both  subject  to  certain  specific 
assumptions  regarding  the  structure  of  the  systems  of  galaxies  and  regarding 
the  distribution  of  peculiar  velocities.  Also,  both  tests  are  optimal  with  regard 
to  the  alternative  hypotheses  that,  if  a  system  of  galaxies  is  unstable,  then  in 
addition  to  random  peculiar  velocities  conforming  with  the  definition  of  stability, 
the  members  of  the  system  have  a  nonrandom  Hubble-type  component,  either 
away  from  the  center  or  towards  it,  whose  magnitude  is  proportional  to  the  dis¬ 
tance  between  the  center  of  the  system  and  the  particular  galaxy.  The  coeffi¬ 
cient  of  proportionality  0  is  analogous  to  the  Hubble  constant. 

The  first  test  criterion  T  is  optimal  on  the  assumption  that,  if  the  systems 
studied  are  unstable,  then  they  are  characterized  by  the  same  value  of  0.  The 
second  test  criterion  Z  has  optimal  properties  with  regard  to  the  assumption 
that  the  systems  of  galaxies  studied  may  have  different  coefficients  of  instability 
0.  Because  the  two  tests  deal  only  with  velocities  and  apparent  magnitudes  of 
the  galaxies  and  are  independent  from  dynamical  considerations,  they  are 
termed  kinematical  tests. 

(iii)  The  basic  assumption  conditioning  the  validity  of  the  kinematic  tests  is 
that,  while  the  objects  for  which  the  data  are  available  may  have  been  purposely 
selected  for  their  brightness  or  position  or  some  other  characteristic,  this  selec¬ 
tion  did  not  affect  the  conditional  distribution  of  the  radial  velocity. 

(iv)  The  application  of  the  two  tests  to  empirical  data  regarding  seven 
clusters  and  fifteen  groups  indicated  several  cases  where  the  hypothesis  of  sta¬ 
bility  can  be  rejected  with  little  risk  of  error.  Unfortunately,  conclusions  of  this 
kind  depend  heavily  on  the  basic  assumption  (iii)  and,  more  specifically,  on  the 
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objects  that  one  decides  to  consider  members  of  a  given  system  and  on  the 
systems  that  one  is  inclined  to  consider  as  satisfying  the  hypothesis  (iii). 

(v)  An  outstanding  example  of  this  kind  is  the  Coma  Cluster.  If  one  accepts 
that  all  the  50  objects  presumed  to  be  members  of  the  Coma  Cluster  are  really 
its  members,  then  there  is  little  doubt  that  the  cluster  is  expanding.  On  the 
other  hand,  if  one  excludes  the  “outliers”  there  appears  to  be  some  slight  evi¬ 
dence  of  contraction.  It  is  possible  that  the  “outliers”  are  located  in  the  outer 
regions  of  the  cluster,  which  expands,  while  the  inner  part  of  the  cluster  con¬ 
tracts. 

Similarly,  treating  jointly  all  the  15  groups  of  galaxies  one  sees  no  evidence  of 
instability.  On  the  other  hand,  if  one  removes  from  the  analysis  five  groups  for 
which  there  is  some  evidence  of  selection  for  similarity  of  redshift,  there  is  a 
strong  indication  of  rapid  expansion. 

(vi)  The  basic  numerical  work  was  conducted  on  the  a  priori  assumption 
that  the  systems  of  galaxies  are  essentially  free  of  absorbing  matter.  A  study  of 
the  asymptotic  power  of  the  tests  indicated  that  in  the  absence  of  absorption 
the  tests  would  have  little  chance  of  detecting  instability  even  if  this  were  char¬ 
acterized  by  large  values  of  0.  Since  some  indications  of  instability  were  obtained, 
this  suggests  not  only  that  these  systems  are  unstable  but  also  that  they  contain 
absorbing  matter.  The  suggestion  that  absorbing  matter  is  present  within 
cluster  and  groups  of  galaxies  is  reinforced  by  the  positive  correlation  between 
the  absolute  values  of  the  estimates  of  0  and  the  mean  radial  velocities  of  par¬ 
ticular  systems.  This  correlation  is  exhibited  in  figure  4.  In  the  first  panel  of  this 
figure  each  point  corresponds  to  the  actually  computed  p.  and  |0|  for  each  of  the 
seven  clusters.  In  the  second  panel  the  exhibited  crosses  are  “normal  points” 
corresponding  to  groups  classified  in  intervals  of  1000  km/sec  of  their  radial 
velocities.  Here  each  ordinate  represents  the  weighted  mean  of  the  particular 
values  of  |0,|.  The  presence  of  a  positive  correlation  is  unmistakable,  which  leads 
us  to  the  final  conclusion. 

(vii)  With  reference  to  equation  (60)  it  is  seen  that,  if  either  0  or  r  were  equal 
to  zero,  then  there  would  be  no  correlation  between  0  and  p..  Since  r  cannot  be 
negative,  positive  correlation  is  possible  only  when  0  >  0  and  negative  correla¬ 
tion  only  when  0  <  0.  The  values  of  0i  obtained,  some  positive  and  some  nega¬ 
tive,  suggest  that  some  of  the  systems  are  expanding  and  some  contracting  so 
that  0  might  be  considered  as  a  random  variable  with  some  fixed  mean,  probably 
close  to  zero,  and  a  fixed  variance.  On  this  assumption,  there  would  be  observed 
changes  in  signs  of  the  estimates  0,-,  an  increase  in  the  conditional  dispersion, 
proportional  to  (1  +  r  D/a),  and  a  positive  correlation  between  |0»|  and  p., 
everything  exactly  as  in  the  empirical  results  and  as  reflected  in  figure  4.  This 
particular  hypothesis,  of  positive  r  and  a  random  0,  appears  to  be  worthy  of 
closer  consideration.  In  studying  it,  the  existence  or  nonexistence  of  a  correlation 
between  the  color  index  of  the  galaxies  and  their  radial  velocities  is  likely  to  be 
decisive.  For  expanding  systems  this  correlation  should  be  positive  and  for 
contracting  negative. 
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Correlation  between  the  estimated  absolute  value  of  the  expansion 
coefficient  0  and  the  mean  radial  velocity  of  the  system. 

The  vertical  bar  is  of  length  2 <rg. 

The  upper  panel  refers  to  clusters,  the  lower  panel 
refers  to  groups  (normal  points). 
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